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I  .   INTRODUCTION 

This  paper  outlines  in  detail  the  work  achieved  to  date 
at  the  Naval  Postgraduate  School  in  the  area  of  underwater 
shock  loading  of  plates.   The  objective  of  this  work  is  three 
fold:   first  to  obtain  experimental  data  on  the  dynamic 
response  of  simple  structures  subject  to  underwater  shock, 
secondly  to  compare  the  experimental  results  with  the 
predicted  results  based  on  the  existing  computer  code,  and 
thirdly  for  the  school  as  a  whole  to  acquire  the  knowledge 
and  technical  expertise  necessary  to  conduct  these  types  of 
experiments.   An  anticipated  result  of  the  latter  is  to 
identify  the  problem  areas  and  base  technology  needs  in 
experiment  and  prediction  methods.   This  paper  discusses  the 
first  underwater  shock  test,  which  was  an  investigation  into 
the  associated  phenomena  of  gross  plate  response  and  the 
tripping  of  rectangular  stiffeners. 

The  tripping  effect  is  a  lateral,  torsional  instability 
of  the  stiffener  as  it  becomes  suddenly  unstable  and  fails 
under  a  load.   Tripping  is  also  viewed  as  a  buckling  and 
warping  of  the  stiffener.   In  either  case,  the  response  of 
such  structures  as  a  stiffened  flat  plate  or  a  cylinder  with 
ring  stiffeners  will  change  dramatically  when  the  functional 
character  of  the  stiffeners  is  reduced  after  they  have 
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tripped.   The  need  for  this  study  and  its  applicability  to 
the  Navy  is  self-evident  as  stiffeners  are  incorporated  in 
the  structural  design  of  all  ships  and  submarines.   The 
results  of  a  previous  series  of  tests  on  ship-type  grillages 
have  been  described  as  "clearly  demonstrating  the  signi- 
ficance of  lateral- torsional  instability  (tripping)  as  a 
primary  ductile  failure  mode  for  ship  structure"  [Ref.  1]. 
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II .   DESIGN  AND  FABRICATION  OF  THE  TEST  STRUCTURE 

The  purpose  of  this  work  was  to  conduct  a  series  of 
experimental  tests  where  plates  of  various  stiffener  types 
would  be  loaded  by  an  underwater  shock  wave.   The  plate  was 
geometrically  similar  to  a  ship's  stiffened  hull--or 
grillage — and  air  backed.   This  required  a  backing  structure 
that  was  watertight  and  strong  enough  to  be  used  throughout 
the  test  series.   The  plate  was  subjected  to  large  deforma- 
tions, well  into  plastic  regime,  to  ensure  a  good  tripping 
effect  in  the  stiffeners.   As  an  aid  in  setting  up  and 
validating  each  test,  it  was  necessary  to  model  the  test 
plate  using  a  finite  element  computer  code.   Therefore,  to 
ease  the  modeling  requirements,  the  structure  was  kept 
simple  throughout  the  design  process.   Well  defined 
geometric  and  material  properties  and  good  boundary 
conditions  were  specified.   As  in  any  experimental  test,  the 
number  of  uncertainties  were  kept  as  few  as  possible. 

Needed  for  the  plate  material  was  a  conwnon  and  easily 
machinable  metal  with  the  material  properties  of  being 
initially  isotropic  and  exhibiting  very  little  strain 
hardening.   6061-T6  alaminum  was  selected  because  it  has 
these  properties  and  is  readily  available  in  the  sizes 
required.   It  is  easily  machined  and  can  be  welded  as  well. 
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Upon  receipt  of  the  material,  a  section  of  the  blank  was 

removed  and  tension  test  specimens  made  from  it.   The 

uniaxial  tension  tests  were  conducted  to  characterize  the 

6061-T6  aluminum  as  having  a  yield  stress  (a  )  of  43200  psi 

and  an  ultimate  stress  (a  ^  .)     of  44900  psi.   The  recorded 

^  ul  t  '^ 

elongation  is  11.5%.   Data  from  these  tests  are  presented  in 
Appendix  A. 

As  for  the  boundary  conditions  of  the  test  panel,  a 
clamped  arrangement  has  proven  in  past  shock  tests  to  be 
most  effective.   The  integral  edge  arrangement  reported  by 
E.  A.  Witmer  and  R.  Wu  [Ref.  2]  was  adopted.   This  involved 
machining  the  test  panel  out  of  larger  blank  stock  and 
leaving  a  massive  edge  on  all  four  sides.   Previous  designs 
for  a  clamped  boundary  of  a  panel  to  be  shock  tested  had 
employed  such  schemes  as  serrated  clamps,  massive  bolts,  or 
hardened  faces,  all  of  which  had  slipped  to  a  degree.   It 
was  this  final  design  of  an  integral  edge  which  exhibits 
essentially  no  slip  during  loading  [Ref.  3]. 

The  stiffeners  were  likewise  integrally  machined  from 
the  aluminum  blank  to  eliminate  any  uncertainty  that  would 
arise  from  a  welding  procedure.   All  corners  and  inter- 
sections were  rounded  to  remove  any  stress  concentrators. 

To  model  a  typical  ship's  grillage  requires  the  use  of 
scaling  laws.   Based  on  discussions  with  Dr.  R.  P.  Daddazio 
of  Weidlinger  Associates,  Inc.,  there  are  two  important 
parameters,  B  and  X,  for  scaling  of  the  grillage. 
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;b/t)  ,j  a^/E 


and     \    =     (a/kiT)  )   a  /E      (eqn  2.1) 


The  two  parameters  are  defined  such  that  B  is  a  ;T\easure 
of  the  plate's  slenderness  ratio  and  X  is  a  measure  of  the 
slenderness  of  the  longitudinal  stiffeners.   Other  variables 
are  : 

a   =  material  yield  strength 

E   =  Young's  modulus 

b   =  spacing  of  the  longitudinal  stiffeners 

t   =  plate  thickness 

k   =  radius  of  gyration  of  longitudinal  stiffener  acting 

with  an  assumed  effective  width  of  plating  denoted  by 
b 


a   =  spacing  of  the  transverse  frames 

b   =  effective  width  of  plating  given  by  b   =  b  (- 


■) 


[Ref .  1]  . 

The  final  dimensions  of  the  test  panel  are  as  shown  in 
Figure  2.1.   A  calculation  of  the  B  and  X  as  a  check  on  the 
design  is  necessary.   From  reference  4  a  survey  of  typical 
ship's  grillages  yields: 


1  <  B  <  4.  5 


and 


.15  <  X  <  .9 


Approximate  engineering  reference  text  values  for  the 
6061-T6  yield  strength  and  elastic  modulus  were  used  in  this 
design  process  because  the  stock  material  was  not  on  hand 
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Fig.     2.1      Test    Panel 
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for  testing.   For  the  following  calculations,  a^  =  40  ksi 
and  E  =  10  X  10^  psi  . 

For  the  rectangular  stiffener  in  Figure  2.2 


b  \  ^  _  f 6  inches 

"  ^tN  E 


40000psi 


0.1875  inches^  ^  10  xlO   psi 


1/2 


=  2.02 


eqn  2.2 


^e  =  ^(|  -  72^  =  (^  inches)  (  ^702  "  " 7^2 


)  =  4.47  inches 

(eqn  2.3) 


Neutral  Axis  Location, 


(0.5  +  ^'-^^ '^^)  (0.1875)  (1.0 


na        (0.1875)  (4.47  +  1) 


=  0.109  inches    (eqn  2.4 


1^1.^^  =  t4-  (4.47)  (0.1875) ^  +  (4.47)  (0.  1875)  (0.109) ^  (eqn  2.5) 
px  a  tie    X  z 


0.0124  in^ 


Igtiffener  =  IT  (0-1875)  (1)^  +  (0 . 1875 )  ( 1 )  (  0 . 48 5 ) ^    (eqn  2.6 


0.0597  in^ 


Radius    of    Gyration, 


^  =   ^i 


)   I  /0.0124    +    0.0597    \^^'^         r,    nc^     ■       u  ,  0-7 

J  A      =      (-(0-.1875)  (4.47  +  1)  )  =0.265    inches  .(eqn    2.7 
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Figure  2.2   Stiffener  Dimensions 


^=  (^'Je^=  ''  '"='" 


(0. 265 


nes     /  4U  yUUU  psi 
in)  (t7)  ho  X  10^  psi 


1/2 


=  0.9 


(eqn  2.8) 


It  is  noted  that  6=2  and  X  =  .9  fall  within  the  range 
of  a  typical  ship's  grillage  and  the  test  panel  is  an 
acceptable  model . 

The  final  design  is  a  panel  18  inches  by  12  inches  and 
3/16  inch  thick  machined  from  a  2  inch  thick  aluminum  blank 
measuring  33  inches  by  27  inches.   There  are  two  stiffeners 


18 


located  symmetrically  about  the  centerline.   This  is  to 
permit  a  computer  modeling  of  only  one  quarter  of  the  test 
panel  which  is  a  great  savings  in  computational  costs.   The 
stiffeners  are  on  the  exposed  external  surface  so  that  the 
loading  conditions  at  the  center  will  be  compressive  in 
nature;  this  is  a  requirement  for  tripping  to  occur. 

The  test  structure  is  made  in  such  a  way  that  the  plate 
may  be  simply  turned  over  and  secured  to  the  backing 
structure  should  it  be  desired  to  test  internally  stiffened 
plates.   A  bolting  arrangement  between  the  test  panel  and 
the  backing  structure  was  designed  to  aid  the  clamped 
boundary  conditions  and  to  ensure  a  watertight  seal 
throughout  the  loading.   The  final  design  is  a  double  row  of 
64  bolts  all  around.   They  are  1  inch  in  diameter,  A325  high 
strength  structural  steel  in  a  f r iction- type  connection 
under  single  shear  loading.   An  eyebolt  is  located  at  each 
corner  for  ease  in  handling  and  rigging. 

The  design  drawings  for  the  test  plate  and  the  backing 
structure  are  in  Appendix  B  and  are  to  be  referred  to  for 
more  precise  details  and  dimensions.   A  presentation  of  the 
calculations  performed  in  the  design  of  the  bolting 
arrangements  will  be  made  here. 

Figure  2.3  shows  a  bolt  group  located  on  the  massive 
aluminum  collar  surrounding  the  test  section.   Bolt  spacing 
is  3  inches. 
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Figure  2.3   Bolt  Group  Spacing 


Maximum  bolt  shear  =  (a   panel) (t  panel) (bolt  spacing) 

=  (40  Ksi)  (3/16  inch)  (3  inches) 
=  22.5  Kips 


Maximum  moment  at  edge  =  1/4  a  t  b 

=  (1/4) (40  Kips) (3/16  inch)^(3  in; 
=  1.05  Kip- in 

This    moment    produces    a    tension    in    the   bolt: 
p=M/e   =    1.05    Kip-in/3    inches    =    0.3516    Kips. 
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To  check  the  effect  of  this  tension  on  bolt  shear,  the 
American  Institute  of  Steel  Construction  Manaul  [Ref.  5]  was 
used  . 

Specification  1.6.3  on  page  5-28  defines  an  allowable 
shear  reduction  factor  due  to  tension  in  the  bolt. 
Shear  reduction  factor  =  (1-f A/T) : 

where  f  =  average  tensile  stress 
A  =  nominal  body  area 
T  =  minimum  bolt  pretension, 
from  Table  1. 23 . 5 
factor  =  (1-. 3516/51)  =  0.993 
Table  1.5.2.1  lists  the  maximum  shear  stress  allowed,  Fv 
=  17.5  Ksi ,  so  the  allowable  shear  =  (17.5  Ksi)(.993)  = 
17.38  Ksi.   Therefore,  the  effect  of  tension  in  the  bolt  is 
negl ig  ible . 

From  the  AISC  manual  page  4-5,  Table  ID,  the  allowable  shear 
load  for  a  friction  connection  using  one  inch  A325  bolts  is 
(2  bolts  in  line  of  shear)  (13.7  Kips)  =  27.4  Kips 
Check  with  the  maximum  bolt  shear,  22.5  Kips  <  27.4 
Kips,  and  find  that  the  bolting  arrangement  is 
acceptable . 

A  second  set  of  calculations  for  a  bolting 
configuration  where  the  plate  is  flipped  over  also 
proves  satisfactory. 
Maximum  moment  at  the  edge  =  1/4 (a  )  (t)  (b)  (t) 
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=    1/4 (40    Ksi) (3/16) (3) (2) 
=    11.25    Kip-in 
Tension    in    the    bolt,    P=M/e    =    11.25/3  =    3.750    Kips 
Shear    reduction    factor    =     (1-fA/T) 

k    =     (1-. 375/51) 
k    =    0.927,    still    consider    k 
neg 1 ig  ible  . 
The    backing    structure    shown    in    Figure    2.4   was    designed 
using    the   methods    described    in     [Ref.    6].       The    entire 
structure    is   made    of    standard    structural    A36    steel,    1.25 
inches    thick.       An    0-ring    gasket    is    fitted    into    the    channel 
machined    in    the    surface    of    the    flange.       A   hull    penetrator 


Figure    2.4      Backing    Structure    During    Assembly 


22 


type  connector  is  fixed  to  the  bottom  center  to  allow  for 
internal  instrumentation.   A  steel  collar  is  welded  around 
the  exposed  penetrator  to  prevent  damage  during  handling  and 
testing.   The  penetrator  is  a  24-pin  model  made  by  SEACON, 
part  no.  XSM-BCR.   The  mating  connector  used  during  the  test 
is  part  no.  XSM-CCP,  with  a  40  foot  length  of  16  gage  cable 
attached.   The  backing  structure  weighs  approximately  758 
lbs.   The  test  panel  weighs  126  pounds,  and  the  nuts  and 
bolts  weigh  111  lbs.,  for  a  total  weight  of  the  test 
structure  of  1,000  pounds.   This  is  obviously  massive  and 
has  proven  strong  enough  to  withstand  repeated  shock  tests. 
Figure  2.5  shows  the  assembled  structure  on  the  workbench. 


Figure  2.5   Backing  Structure  With  Test  Panel 
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Ill .   TEST  GEOMETRY 

The  inelastic  response  of  the  stiffened  panel  which  this 
test  is  to  examine  requires  that  the  plate  undergo  large 
deflections.   These  are  on  the  order  of  four  plate  thick- 
nesses.  To  achieve  this  result,  calculations  must  be  made 
which  specify  the  appropriate  charge  weight  and  standoff 
distance.   A  treatment  of  the  theory  of  underwater 
explosions  will  not  be  attempted  in  this  paper,  rather  a 
prior  knowlege  is  asumed  concerning  the  characteristics  of 
the  pressure  loading  generated  and  the  many  secondary 
phenomena  associated  with  underwater  explosions. 

The  empirical  formulas  used  to  determine  the  charge 
weight  and  standoff  are  those  derived  by  Robert  H.  Cole  in 
his  authoritative  text,  Underwater  Explosions  [Ref.  7]. 

To  determine  the  pressure  at  a  point  as  a  function  of 
time,  an  equation  is  written  in  a  general  form  as 


P(t)  =  Pmax (e) 


-t/9 


(eqn  3.1) 


where  the  maximum  pressure  is 


Pmax  =  K,  (W-*-/^  /R)  "'■  ps  i 


(eqn  3.2; 
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and  the  exponential  decay  constant  is 


e  =  K2(W-'-'^^)  (W^'^^  /  R)  ^  msec 


eqn  3.3) 


While  the  forms  of  these  equations  have  been  accepted  as 
correct  and  invariant,  the  coefficients  K-,  ,  A,  ,  K„  ,  and  A^ 
are  determined  from  data  taken  from  numerous  tests  and  are 
redefined  as  more  accurate  data  is  obtained.   Due  to  the 
fact  that  the  values  will  differ  with  different  types  of 
explosives,  the  first  step  in  determining  a  test  geometry 
would  be  to  specify  the  explosive.   TNT  was  selected  due  to 
its  frequent  use  in  tests  of  this  sort,  and  an  eight  pound 
charge  was  chosen  due  to  its  availability  in  the  Navy 
requisitioning  system.   The  coefficients  for  TNT  are 


K-^  =  22505 


K2=0.058 


A,  =1.  li 


A2=0. 185 


The  values  of  W  and  R  are  the  charge  weight  in  pounds  and 
the  standoff  distance  in  feet,  respectively. 

Now  that  the  incident  pressure  loading  on  the  plate  is 
known  as  a  function  of  time,  the  theories  of  mechanics  of 
structures  may  be  applied  to  determine  the  resulting 
deformation.   Cole  treats  this  topic  through  the  use  of 
energy  methods  and  draws  a  distinction  between  elastic  and 
plastic  considerations  of  yielding  surfaces.   An  equation  is 
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derived  that  provides  a  value  for  the  final  deflection  of  a 
plate,  Z(t^,)  ,  by  equating  the  plastic  work  to  the  kinetic 
energy  acquired.   The  equation  is  given  below,  but  the 
reader  is  referred  to  Cole  for  its  derivation  and  the 
definition  of  its  variables. 


Z(t^)  = 


/°oS 


2m 

o  h 
o 


^'"    -^JA/5 


^"^  a  E^/'2  (eqn  3.4) 


The  preceding  equations  could  have  been  used  in  an 
iterative  fashion  to  determine  the  required  charge  weight, 
charge  type  and  standoff  to  produce  the  necessary 
deformation  (the  four  plate  thicknesses  is  4  x  0.1875  =  0.75 
inches) .   However,  this  method  of  calculation  was  not 
pursued.   Rather,  a  computer  analysis  was  made  using  a 
finite  element  model  of  the  test  plate.   This  is  a  much  more 
powerful  method  of  analysis  in  terms  of  time  and  accuracy. 
Moreover,  in  addition  to  specifying  plate  deformations, 
regions  of  maximum  stress  and  strain  and  the  stiffener 
response  may  be  investigated,  dependent  upon  the  computer 
code  used.   The  strains  predicted  by  the  computer  output 
were  used  to  set  the  instrumentation  levels  at  the  test 
site.   This  is  one  example  of  the  benefits  to  be  gained  in 
using  such  a  finite  element  model.   A  discussion  of  this 
computer  model  and  the  results  obtained  is  the  topic  of  the 
next  chapter . 
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The  analysis  indicated  the  best  standoff  distance 
between  the  eight  pound  charge  and  the  test  plate  was  nine 
feet . 

As  in  any  experimental  investigation,  the  test 
environment  must  be  as  ideal  as  physically  possible.   The 
test  conditions  desired  were  that  they  be  repeatable,  that 
they  concentrate  on  the  effect  of  the  shock  wave  loading, 
and  that  they  remain  clean  and  free  from  any  of  the 
secondary  effects  associated  with  underwater  explosions.   An 
attempt  was  made  to  eliminate  or  minimize  such  phenomena  as 
bulk  cavitation  and  cavitation  closure,  the  reloadings  from 
the  explosive  gas  bubble  pulse  and  bubble  migration,  and  the 
surface  cutoff  and  bottom  reflection.   To  achieve  all  of 
these  criteria  requires  a  combination  of  correct  test 
geometry  and  limiting  the  elapsed  time  of  the  investigation 
to  a  small  window.   Most  of  these  phenomena  are  considered 
late  time  effects  and  are  eliminated  by  only  looking  at  the 
plate  response  during  the  first  few  milliseconds  of  pressure 
loading.   Furthermore,  this  early  time  restriction  allows 
for  the  modeling  of  the  fluid/structure  interaction  in  a 
most  simple  and  convenient  manner.   The  plate  response  may 
be  determined  by  approximating  the  loading  as  a  plane 
pressure  wave.   To  consider  the  late  time  or  even  inter- 
mediate characteristics  of  the  fluid/structure  interaction 
would  require  entirely  different  approximations  [Ref.  8]. 


27 


One  of  these  late  time  effects  which  can  be  destructive 
in  nature  is  the  reloading  generated  by  the  explosive  gas 
bubble.   There  are  two  distinct  phenomena  caused  by  the  gas 
sphere:   an  expansion  and  contraction  cycle  which  generates 
additional  pressure  pulses  and  thereby  reloads  the 
structure,  and  the  tendency  for  the  bubble  to  migrate 
towards  a  structure  if  it  is  close  enough,  and  to  then 
collapse  upon  it  and  do  destructive  work.   Fortunately, 
empirical  equations  have  been  determined  which  permit  the 
calculation  of  the  bubble  radius  as  a  function  of  time  and 
the  time  of  the  first  closure  pulse.   These  two  equations 
are  similar  in  form  to  the  pressure  equation.   The 
coefficients  are  determined  by  the  type  of  explosive  and  the 
variables  are  the  charge  weight  and  the  charge  depth.   The 
general  equation  and  its  associated        coefficients  for 
TNT  are  as  follows: 


Time  of  first  closure  pulse   T  =  K 


W 


1/3 


5  (D+33)^/^ 


sec   (eqn  3.5) 


Maximum  bubble  radius   A     =  K^ 

max     6 


W 


1/3 


(D+33) 


1/3 


ft 


(eqn  3.6) 


K^  =  4.268 


Kg  =  12.672 
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An  eight  pound  charge  weight  with  a  nine  foot  standoff 
distance  has  already  been  selected  to  produce  the  desired 
deformation.   W  =  8  lb.  and  charge  depth  in  feet  are  run 
iteratively  in  these  equations  to  satisfactorily  meet  two 
conditions.   First,  the  time  of  the  closure  pulse  is  to  be 
relatively  late,  well  beyond  the  first  few  milliseconds. 
Secondly,  the  charge  depth  needs  to  be  such  that,  as  the 
bubble  expands  to  its  maximum  radius,  it  will  break  the 
surface  of  the  water  and  vent  to  the  atmosphere  before  it 
contacts  the  plate.   This  venting  action  is  considered 
instantaneous  and  the  bubble  is  prevented  from  expanding  or 
contracting  further  if  it  breaks  the  surface  when  it  is  near 
its  maximum  radius.  A   ratio  of  charge  depth  to  maximum 
radius  in  the  neighborhood  of  .50  to  .75  will  ensure  a  good 
venting  action.   A  charge  depth  of  four  feet  was  determined 
to  produce  the  desired  results.   The  calculations  are  shown 
below. 


T  =  (4.268) 


8 


1/3 


(4+33) 


5/6 


=  0.421  sec 


^max  =  (12.67) 


8 


1/3 


(4+33) 


1/3 


=  7.60  ft 


4A/7.60  ft  =  0.526 
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A  charge  depth  of  four  feet  is  53%  of  maximum  bubble 
radius  and  it  will  therefore  vent. 

Another  underwater  explosion  phenomenon  which  can  be 
eliminated  through  the  proper  test  geometry  is  the  bulk 
cavitation  effect.   The  region  of  water  that  will  cavitate  is 
a  function  of  depth  and  range  from  the  explosion  source.   The 
cavitation  is  created  when  the  reflected  wave  from  the 
surface  passes  through  the  water  directly  behind  the  primary 
wave  front.   Figure  3.1  depicts  the  generation  of  bulk 
cavitation  (courtesy  of  Weidlinger  Associates) . 

The  effect  this  cavitated  region  has  on  the  structure  is 
that  the  plate  will  experience  essentially  no  pressure 
loading  while  it  is  surrounded  by  cavitated  water.   However, 
when  the  combined  forces  of  atmospheric  pressure  and  the 
weight  of  the  water  above  the  region  overcome  the  cavitating 
forces,  the  cavitation  suddenly  closes  back  up  to  generate  a 
reloading  on  any  structure  within  the  region.   Although  the 
pressure  associated  with  this  cavitation  closure  can  be 
calculated  and  has  been  successfully  done  by  Weidlinger 
Associates  in  a  finite  element  code  called  CUE,  these 
calculations  are  extremely  involved.   The  best  solution  for 
the  experimental  test  in  question  would  be  to  remove  the  test 
structure  from  the  cavitated  region  entirely.   This 
necessitates  the  use  of  equations  which  define  the  extent  to 
which  cavitation  will  occur,  based  on  charge  size,  type,  and 
location  . 
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Figure  3.1   Development  of  Cavitation  as  Described  by  CUE 

The  equations  ncessary  to  calculate  this  region  of  bulk 
cavitation  are  substantial  in  size.   Computation  of  this 
region  by  computer  is  advisable  and  was  done  by  means  of 
another  program  written  by  Weidlinger  Associates,  Inc.  for 
Tektronix  4051/2.   The  graphic  output  from  this  program  is  a 
range  versus  depth  profile  of  the  cavitated  region  and  is 
shown  in  Figure  3.2. 

The  result  of  all  calculations  is  a  vertical  charge/plate 
orientation  with  the  charge  located  directly  over  a  flat, 
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Figure    3.2      Region    of    Bulk    Cavitation 
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horizontal  plate.   The  eight  pound  TNT  charge  would  be 
suspended  at  a  depth  of  four  feet  and  the  test  structure  nine 
feet  beneath  it  at  a  depth  of  thirteen  feet.   Theoretically, 
this  test  geometry  would  create  a  test  environment  which  is 
clean  from  those  effects  which  are  burdensome  to  calculate 
and  reproduce  while  still  producing  the  desired  deformations 
in  the  plate. 

The  previously  mentioned  computer  program  written  for  the 
Tektronics  4051/2  by  Weidlinger  Associates  is  part  of  a 
library  of  programs  available  at  the  Naval  Postgraduate 
School.   A  second  program  on  this  tape  library  is  called 
"Undex  Parameters"  and  outputs  numerical  data  for  most  of  the 
phenomena  of  interest  in  an  underwater  explosion.   Copies  of 
the  output  for  this  specific  test  geometry  and  the  pressure 
history  graph  generated  are  presented  on  the  following  pages. 
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UHDEX   PARAMETERS   
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TEST  GAGE 

1 

CHARGE  WEIGHT 

3.8008 

LBS.  TNT 

STANDOFF 

9.8888 

FT 

CHARGE  DEPTH 

4.8888 

FT 

MODEL  DEPTH 

13.8888 

FT 

HORIZONTAL  RANGE 

3.3088 

FT 

PEAK  PRESSURE 

3947.5134 

PS  I 

DECAY  CONSTANT 

0.1573 

MSEC 

SHOCK  FACTOR 

0.3143 

SHOCK  FACTOR  WITH  CUTOFF 

8.3143 

KEEL  SHOCK  FACTOR 

8.3143 

KEEL  SHOCK  FACTOR  WITH  CUTOFF 

3.3143 

ENERGY 

226.9166 

PSI-IH 

ENERGY  WITH  CUTOFF 

226.9166 

PS  I -IN 

IMPULSE 

0.7656 

P3I-SEC 

IMPULSE  WITH  CUTOFF 

8.7656 

PSI-SEC 

CUTOFF  TIME 

1.6868 

MSEC 

BUBBLE  PERIOD 

0.4302 

SEC 

BUBBLE  RADIUS 

7.5625 

FT 

BUBBLE  MIGRATION 

7.2622 

FT 

PRESSURE  REDUCTION  FACTOR 

8.1816 

BUBBLE  PULSE  PRESSURE 

689.7383 

PSI 

BULK  CAUITATION  DEPTH 

1.5758 

FT 

TIME  OF  CAUITATION  CLOSURE 

172.5546 

MSEC 

(CUTOFF  TIME)/<DECAY  CONSTANT) 

10.1789 

Figure    3.3      UNDEX    Parameters    for    Test    Geometry 
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IV.   THE  EPSA  MODEL 

The  finite  element  computer  code  used  to  predict  the 
gross  plate  reponse  was  the  Weidlinger  authored  Elasto- 
Plastic  Shell  Analysis  (EPSA)  code.  The  theory  underlying 
the  analytic  expressions  used  by  the  code  may  be  found  in 
the  EPSA  User's  Manual  [Ref.  12].   The  version  installed  at 
NPS  is  run  on  a  VAX-11/780  machine  using  the  EPSA  Guide 
[Ref.  9].   An  interface  between  EPSA  and  the  PATRAN-G 
graphics  software  was  the  subject  of  recent  thesis  work  at 
NPS  [Ref.  10].   This  eased  the  task  of  interpreting  the 
numerical  output  and  has  created  a  powerful  tool  for  the 
analysis  of  underwater  shock  loading  of  structures  (both 
cylinders  and  flat  plates) .   Appendix  C  presents  the 
necessary  commands  and  procedures  for  executing  an  EPSA  job 
on  the  NPS  installation.   Also  presented  here  are  examples 
of  file  editing  and  the  use  of  library  files. 

EPSA  is  a  central  difference,  finite  element  numerical 
scheme  which  will  output  nodal  point  and  element  response 
quantities  at  the  specified  number  of  time  steps. 
Computation  starts  when  the  incident,  spherical  shaped 
pressure  wave  loads  the  first  node  it  reaches  based  on  the 
input  geometry.   Output  quantities  are  the  incident 
pressures,  incident  velocities,  and  the  total  pressure  at 
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the  requested  nodes.   Also  output  are  the  requested  element 
stresses  and  strains,  nodal  velocities  and  displacements. 
It  was  these  displacements  which  were  used  to  check  the  test 
geometry  to  ensure  the  necessary  deformation  was  occurring 
during  the  initial  loading.   Also,  as  part  of  the  other  pre- 
shot  calculations,  the  strains  were  output  at  the  locations 
of  the  strain  gages  on  the  plate.   The  maximum  value  of 
strain  was  used  to  set  the  instrumentation  levels  prior  to 
the  test . 

Because  this  test  assumes  that  theory  can  only 
accurately  predict  the  plate  response  for  the  first  few 
milliseconds  of  shock  wave  loading,  EPSA  was  run  for  this 
limited  time  frame.   Modifications  were  made  to  the  versions 
of  EPSA  at  MPS  so  that  only  the  early  time  effects  are 
considered.   The  fluid-structure  interaction  used  only  a 
plane  wave  approximation;  the  later  time  effects  within  the 
code  (the  virtual  mass  approximation  in  the  doubly 
asymptotic  approximation)  are  turned  off.   Additionally, 
since  EPSA  was  originally  written  to  model  cylindrical 
shells,  either  with  or  without  internal  stiffeners, 
modifications  were  made  to  permit  the  modeling  of  flat 
plates  with  external  stiffeners. 

The  procedure  for  creating  an  EPSA  input  deck  is  covered 
in  the  User's  Manual.   Only  those  input  calls  which  are 
different  will  be  discussed.   Figure  4.1  is  the  input  deck 
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Figure  4.1   EPSA  Input  Deck  for  Preshot  Analysis 
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used  for  the  present  analysis.   Note  that  the  timestep  size 
is  5  X  10   seconds  and  the  program  is  to  run  for  4,000  time 
steps.   Therefore,  the  real  time  loading  of  the  plate  which 
is  output  is  2  milliseconds,  and  a  run  of  this  time  size 
requires  1  hour  and  32  minutes  of  CPU  time  on  a  VAX  11/780 
computer . 

The  discretization  of  the  plate  and  the  axes  for  the  i-j 
counting  scheme  discussed  in  the  User's  Manual  are  shown  in 
Figure  4.2.   A  quarter  model  of  the  plate  was  generated  to 
conserve  computer  resources.   The  quarter  plate  is  16 
elements  by  10  elements  with  187  nodes;  each  element  would 
physically  measure  0.52  inches  by  0.60  inches.   An 
application  of  the  Courant  stability  criterion  for  finite 
element  codes  shows  the  timestep  size  to  be  sufficiently 
small  to  avoid  numerical  instability.   This  is  also  shown  in 
the  Figure  4.2. 

When  modeling  a  structure  and  using  the  EPSA/PARTRAN-G 
interface,  care  must  be  taken  so  that  the  origin  of  the 
sheet  is  such  that  side  four  is  either  clamped  or  a  rigid 
end  cap  of  a  cylinder.   This  is  due  to  the  method  by  which 
deformed  meshes  are  generated  by  subtracting  the  rigid  body 
motion  of  the  structure  from  total  displacements.   The 
subroutine  NEUDISP  in  the  interface  module  subtracts  the 
motion  of  the  nodes  on  side  four  from  the  other  nodes  along 
corresponding  rows.   This  anomaly  only  affects  the  PATRAN 
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deformed  mesh;  the  EPSA  results  will  be  correct  regardless 
of  origin  location. 

To  execute  the  plane  wave  approximation,  the  following 
variables  for  the  input  of  Section  VI  of  the  User's  Manual 
must  be  set  to  the  values  listed: 
NCRTOT  =  -1 
MAXTOT  =  1 
IFLSYM  =  1 

This  satisfies  GO  TO  loops  which  pass  over  the  sections 
of  the  code  which  calculate  the  virtual  mass  approximation 
in  the  doubly  asymptotic  approximation.   Additionally,  to 
model  a  flat  plate,  the  input  value  for  the  variable  CURV(2) 
found  in  Section  IV  must  be  a  small  non-zero  number, 
approximately  1.0  x  10 

The  modifications  of  EPSA  which  permit  it  to  treat  the 
case  of  external,  rather  than  internal,  stiffeners  are 
contained  within  the  code.   Each  time  the  other  stiffener 
type  is  desired,  changes  must  be  made  to  the  code  itself. 
This  involves  an  editing  of  the  Fortran  version  of  EPSA, 
followed  by  a  compiling  and  a  linking  to  a  library  file. 
Examples  of  this  procedure  are  found  in  Appendix  C.   The 
EPSA  User's  Guide  also  outlines  this  procedure  [Ref .  12] . 
The  change  required  is  in  the  subroutine  RDSTIF.   Twenty- 
four  lines  from  the  beginning,  change  CX=CX  to  CX=  -CX  when 
it  is  desired  to  change  from  an  internal  stiffener  to  an 
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external  stiffener.   CX  is  the  distance  from  the  centerline 
of  the  plate  cross-section  to  the  stiffener  centroid.   The 
type  of  stiffener  being  modeled  by  EPSA  can  be  determined 
from  the  output  file.   In  the  first  section  following  the 
input  readback  is  a  section  listing  sheet  stiffener  data.   A 
negative  value  for  CX  confirms  the  desired  case  of  an 
external  stiffener. 

Another  update  to  the  EPSA  at  NPS  which  is  not  reflected 
in  the  User's  Manual  permits  the  option  to  output  strains  at 
the  top  of  a  rectangular  stiffener.   To  call  for  this 
output,  the  variable  ITYPE  found  in  Section  XIII  is  set  to 
11.   A  copy  of  the  changes  made  within  the  code  to  the 
subroutine  CPTSTRN  may  be  found  on  page  94  of  reference  12. 

An  additional  modification  to  EPSA  allows  for  the  input 
of  a  discrete  pressure  history  for  the  loading.   The 
discrete  form  of  pressure  field  representation,  as  opposed 
to  an  exponentially  decaying  curve,  is  selected  by  setting 
the  variable  ICHRG  =  2  in  Section  XVI.   Card  3  now  reads 
NSHAPE  (up  to  a  maximum  of  J  =  25  points) .   Card  4  reads 
pairs  of  TC(J)  and  PC(J),  where  TC  is  an  array  of  times  in 
milliseconds  after  the  arrival  of  the  incident  pressure  wave 
and  PC  is  the  array  of  pressure  which  corresponds  to  the 
times  in  array  TC .   The  final  two  values  on  card  4  are 
THEXP2  and  TCUT2.   An  exponentially  decaying  tail  may  be 
added  to  the  discrete  pressure  history.   THEXP2  is  the  time 
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constant  of  exponential  decay  for  this  tail  (corresponds  to 
a  theta  value  in  milliseconds) .   If  a  decaying  exponential 
tail  is  not  desired,  set  THEXP2  <_  0.0.   Figure  4.3  is  an 
example  of  this  option.   TCUT2  is  the  value  of  surface 
cutoff  by  the  negative  reflected  wave  from  the  surface  with 
time  in  milliseconds. 

One  other  unique  feature  about  executing  an  EPSA  flat 
plate  model  is  that  there  is  no  need  to  create  a  virtual 
mass  array  file.   This  fact  is  covered  in  the  User's  Manual 
by  the  example  of  the  statically  loaded  beam.   When  modeling 
a  cylinder  and  employing  the  VMA,  a  virtual  mass  array  must 
be  created  in  machine  language  by  using  either  the 
interactive  program  ACESION,  or  the  batch  BASESION.   For  the 
EPM  models  (explosive  power  meter),  this  file  is  usually 
named  EPM.VMA.   When  running  a  flat  plate,  simply  return 
empty  the  call  for  a  virtual  mass  file  name. 

As  mentioned  earlier,  a  discussion  on  the  merging  of 
PATRAN  and  EPSA  and  the  commands  necessary  for  its  use  was 
the  topic  of  a  previous  thesis.   Therefore,  it  will  not  be 
presented  here,  but  a  brief  description  of  it  and  example  of 
its  use  may  be  found  in  Appendix  C.   This  merger  causes  EPSA 
to  create  special  result  files  of  numerical  data  during  eacn 
run.   These  files  are  formatted  such  that  they  may  be  read 
by  PATRAN.   The  neutral  results  file  is  written  on 
FOR019.DAT  and  contains  the  finite  element  mesh  generated  by 
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Figure  4.3   Discretized  Pressure  History  Input 

the  EPSA  input.   The  mesh  is  shown  in  Figure  4.4.   This  view 
assures  the  user  that  the  input  file  is  correct  and  the  mesh 
is  generated  as  desired.   There  are  187  nodes  and  160 
elements  for  the  flat  plate  which  models  the  upper  left 
quarter  of  the  test  panel.   The  elements  are  numbered  in 
Figure  4.4.   Files  FOR016.DAT  and  FOR018.DAT  are  the 
elemental  results  and  nodal  results  files,  respectively,  and 
are  used  for  post-processing.   IN  EPSA,  an  element  result  is 
either  a  stress  or  a  strain  and  a  node  result  is  either  a 
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Figure  4.4   EPSA-Generated  Finite  Element  Mesh 

velocity  or  a  displacement.   The  PATRAN  graphics  used  for 
the  preshot  analysis  were  deformed  mesh  views  of  the  plate 
and  deformation  contour  level  views.   These  were  generated 
by  running  EPSA  for  2  milliseconds  of  real  time  loading  and 
specifying  an  eight  pound  TNT  charge  at  a  nine  foot 
standoff.   These  views  are  presented  in  Figure  4.5  and 
Figure  4.6.   A  view  of  the  stress  contour  levels  is  not 
shown  here  because  the  Von  Mises  stresses  which  were  to  be 
output  by  PATRAN  (EPSA  provides  only  a  ,  a  ,  and  <J„^)  were 
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Figure  4.5   Deformed  Mesh 

not  being  calculated  correctly  and  it  was  not  possible  to 
fix  this  until  after  the  shot.   The  Von  Mises  stresses  are 
shown  in  the  results  section. 

As  pictured  above,  the  deformed  mesh  provides  an 
excellent  analysis  tool;  the  areas  of  greatest  deformation 
are  easily  identified. 

During  the  post-processing,  PATRAN  will  specify  the 
greatest  displacement  and  the  node  where  it  occurs.   For  the 
deformed  mesh  shown  it  is  -.75  inches  at  node  17.   A  word  of 
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Figure  4.6   Deformation  Contours 

caution  to  the  user:   the  deformed  mesh  is  not  scaled,  but 
rather  exaggerated  so  that  the  results  may  be  prominent. 
This  may  be  deceiving  at  times.   Remember,  the  deformed 
geometry  mesh  is  not  proportioned  exactly. 

The  contour  level  plots  are  another  means  of  displaying 
the  displacements  of  the  plate.   Although  they  may  not 
provide  as  clear  an  image  of  the  results  as  the  deformed 
mesh,  they  do  provide  a  more  accurate  analysis  tool.   With 
each  contour  line  is  an  identifying  letter  which  is 
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transcribed  by  the  legend  so  that  the  precise  amount  of 
displacement  may  be  read  for  any  location  on  the  plate  from 
this  one  view.   Here  it  can  be  seen  that  the  area  around 
node  17  is  experiencing  a  displacement  of  at  least  -.739 
inches . 

These  PATRAN  views  were  checked  against  the  EPSA  output 
file  and  the  displacements  were  verified.   The  requisite 
amount  of  plate  deformation  is  occurring  within  the  first 
few  milliseconds  of  plate  loading  (again  this  is  considered 
the  only  time  frame  which  theory  can  predict  acceptably 
well).   From  the  EPSA  printout  of  displacements,  the 
greatest  displacement  during  this  time  is  -.974  inches  at 
1.5  milliseconds. 

The  final  use  of  EPSA  in  the  preshot  analysis  was  the 
identification  of  maximums  strains  for  setting  instrumen- 
tation levels.   To  better  understand  how  this  is  done,  the 
EPSA  output  for  strains  corresponding  to  gage  locations  on 
the  plate  are  shown  in  Figures  4.7,  4.8,  and  4.9.   The 
maximum  strains  are  circled  and  a  summary  of  these  maximum 
strains  at  each  of  the  gage  locations  is  shown  in  Table  I. 
The  strain  gage  locations  are  shown  in  Chapter  V. 
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TABLE  I 
Summary  of  EPSA  Predicted  Strains 


Gage 

EPSA 

Max  Strain 

Time 

Percent 

Number 

Number 
16 

( p^s  train) 
+  7070 

(sec) 
0.80 

(%) 

1 

0.7 

2 

19 

-9530 

1.95 

0.9 

3 

18 

-7210 

2.00 

0.7 

4 

17 

+  4680 

1.75 

0.5 

5 

15 

-3170 

2.00 

0.3 

6 

12 

-10300 

1.90 

1.0 

7 

13 

+16900 

2.00 

1.7 

8 

14 

-10300 

1.80 

1.0 

9 

11 

-55100 

1.  50 

5.5 

10 

10 

-46900 

0.85 

4.7 

11 

8 

-33300 

1.00 

3.3 

12 

9 

-63000 

1.00 

6.3 

13 

6 

+26700 

1.00 

2.7 

14 

3 

-13400 

1.35 

1.3 

15 

4 

+16900 

2.00 

1.7 

16 

7 

* 

• 

* 

17 

1 

-3170 

2.00 

0.3 

18 

2 

-5570 

1.65 

0.6 

19 

5 

-10300 

1.90 

1.0 

20 

20 

+48800 

0.80 

4.9 

21 

21 

-58900 

1.05 

5.9 
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V.   TESTING  AT  WEST  COAST  SHOCK  FACILITY 

With  the  design  and  manufacture  of  the  test  structure 
completed,  a  satisfactory  test  geometry  determined,  and  a 
pre-shot  analysis  accomplished  with  a  computer  model,  the 
next  step  in  the  experiment  was  to  determine  what 
instrumentation  would  be  necessary,  and  where  and  how  to 
attach  it  . 

The  usual  suite  of  instrumentation  for  an  underwater 
shock  test  consists  of  pressure  gages  in  the  free-field  to 
sense  the  pressure  loading  generated  by  the  charge,  strain 
gages  on  the  test  surface  to  output  strain  and  deformation, 
and  velocimeters  and  accelerometer s  on  the  structure  to 
measure  the  velocity  and  acceleration  imparted  by  the 
impulsive  load.   The  decision  on  the  instrumentation  was 
made  after  liaison  with  the  test  facility  because  it  was  a 
function  of  what  they  could  support  with  their  equipment. 

The  West  Coast  Shock  Facility  (WCSF)  was  selected  as  the 
site  for  the  test  due  to  several  factors.   First  and 
foremost  was  that  they  are  the  only  licensed  facility  on  the 
West  coast  where  a  test  of  this  type  can  be  performed  in  a 
controlled  manner.   Moreover,  the  shock  facility  is  located 
only  120  miles  north  of  Monterey  in  the  South  Bay  of  San 
Francisco.   WCSF  is  a  Navy  activity  which  reports  to  the 
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Supervisor  for  Shipbuilding,  Conversion  and  Repair,  San 
Francisco.   The  facility  has  experience  in  conducting  tests 
on  this  small  scale,  although  the  majority  of  their  work  is 
in  the  area  of  shock-qualification  tests  for  the  Navy.   The 
latter  involves  the  use  of  the  standard  Floating  Shock 
Platforms  (FSP)  and  much  larger  charges.   Design  engineers 
and  an  instrumentation  technician  were  available  to  the 
Naval  Postgraduate  School  to  assist  with  the  details  of  a 
viable  solution  to  the  necessary  test  geometry  and  the 
instrumentation  of  the  structure.   Presented  on  the  next 
page  is  a  sketch  of  the  actual  floation  and  rigging  network 
as  conceived  by  the  design  engineers. 

Figure  5.1  shows  the  actual  test  geometry.   Floation  for 
the  test  structure  was  provided  by  the  two  pneumatic  fenders 
(cylindrical  shaped) .   The  strong  back  I-beam  was  removed 
prior  to  the  shot,  as  it  was  used  solely  for  positioning  the 
rig  into  the  water.   A  wire  rope  was  run  from  the  pier  and 
shackled  to  the  backing  structure  as  a  safety  measure.   This 
would  be  used  to  retrieve  the  structure  from  the  bottom 
should  the  explosion  part  the  supporting  lines.   Once  the 
test  rig  was  submerged  and  the  strong  back  removed,  the 
outhaul  was  taken  across  the  water  by  boat  to  a  parallel 
pier  and  drawn  tight  by  a  winch.   The  test  was  conducted 
between  these  two  piers,  approximately  100  feet  from  the 
nearest  one,  and  in  about  50  feet  of  water.   The  tide  was 
high  and  little  to  no  current  was  present. 
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The  shock  facility  could  support  the  test  with  2  FM  tape 
recorders  which  would  allow  for  24  tracks  of  data  recording 
(14  tracks  per  recorder,  12  tracks  for  data,  with  one  track 
for  voice/countdown  and  one  track  for  a  one  KHz  time 
signal) .   It  was  decided  that  the  24  tracks  be  divided  among 
two  pressure  gages,  one  accelerometer ,  and  twenty  one  strain 
gages.   No  velocimeter  was  used  during  the  test  because 
there  was  not  enough  space  inside  the  air-backing  structure 
where  it  would  be  mounted. 

The  accelerometer  used  was  an  Endevco  pei zoresi stive 
gage,  model  2262-2000,  with  a  range  of  +2000  g.   Ordinarily, 
an  accelerometer  would  be  mounted  into  the  back  of  the  plate 
so  the  acceleration  imparted  to  it  may  be  recorded. 
Additionally,  this  acceleration  could  be  integrated  to 
obtain  a  plate  velocity,  and  integrated  again  for  displace- 
ment.  A  calculation  of  the  theoretical  plate  acceleration 
was  made  using  the  Taylor  Flat  Plate  Theory  to  see  if  it 
was  within  the  limit  of  the  gage.   The  procedure  found  in 
[Ref.  13]  was  followed. 
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=  0.01142  ^^^. ^^^    X  12  in/ft  =  0.137  slugs 
in  ^ 


m  per  sq.  ft. 


=  0.137  slugs/1.5  ft^  =  0.0913  — -^i|H 


1/3 


P'o  =  22505  (  ^ 


)  ■^•■^^  =  3814.9  psi 


1/3 
=  (0.058)  (8^'^^)  (  ^-H-  )  "°-^^^  =  0.1532  msec 


pcQ    ^     (2)  (5000  ft/sec)  (0  .1532x10  ^    sec) 
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(0.0913  Ibf-sec^/ft^) 


=  16.77 
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,  ln3  .     ,_  TCTo  io-3s/  In  16.77  . 
(  gTI  )  =  (0.1532x10   )(   16. 77-1  ^ 


=  0.0274x10    sec 
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2(3814.9x144)  (0.1532x10   )  I     0.1532 
(0.137)  (1-16.77)       I  ® 


0.0274) 


-0  .0274 
0. 1532 


=  61.25  ft/sec 


V 


average  acceleration,  a  = 


max  _   16.25  ft/sec 


2.24x10-"  ft/sec 


max    0.0274x10    sec 


or  69,000  g' 

This  is  in  excess  of  the  Endevco  accelerometer ' s  rating 
or  of  any  accelerometer.   The  plate  is  too  light  and  the 
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impulse  too  severe  to  attach  an  accelerometer  to  it.   The 
other  alternative  was  chosen  and  the  accelerometer  was 
screw-mounted  into  the  base  of  the  backing  structure.   In 
this  manner,  responses  of  the  entire  test  structure  were 
monitored,  and  as  discussed  later,  with  good  results. 

Locations  for  the  twenty-one  strain  gages  were 
determined  based  on  the  following  considerations:   symmetry, 
maximum  plate  response,  stiffener  response,  and  finite 
element  discretization.   Gages  were  located  in  the  X  and  Y 
directions  so  that  the  recorded  output  could  be  compared  to 
the  EPSA  output.   Pairs  of  gages  were  positioned  in  the 
center  cavity  of  the  panel  to  record  maximum  plate  response 
and  a  pair  was  positioned  on  the  underside  of  the  plate, 
internal  to  the  backing  structure,  to  check  outer  fiber 
strains  against  an  exterior  pair  of  gages.   The  balance  of 
the  strain  gages  were  divided  among  the  stiffeners  and  the 
plate.   They  were  mounted  on  the  top  of  the  stiffeners  and 
on  the  sides,  both  in  the  Y  direction  only. 

Figure  5.2  shows  the  strain  gages  fixed  to  the  exterior 
surface  of  the  test  panel.   The  strain  gages  used  were  Micro 
Measurements  type  EA-06-250BG-120 ,  120  ohm,  rated  at  3%  to 
5%  strain  limit  in  tension  or  compression.   The  cement  used 
was  Permabond  910.   For  a  complete  description  of  the 
bonding  technique,  refer  to  Appendix  D. 


Figure  5.2   Photograph  of  Strain  Gages  Bonded  to  Test  Panel 

Figure  5.3  on  the  following  page  is  another  view  of  the 
strain  gage  locations  and  the  position  of  the  accelerometer . 

An  adhesive  sealant  was  required  to  secure  the  strain 
gage  lead  wires  to  the  plate.   On  this  aspect,  the  advice  of 
the  engineers  at  the  Underwater  Explosion  Research  Division 
of  DTNSRDC  Norfolk,  VA  was  followed.   A  saltwater  resistant 
epoxy  was  applied  directly  on  the  gages  and  the  wires.   This 
smoothed  out  while  setting  to  allow  for  good  hydrodynamic 
flow.   The  material  is  a  dichromate  sealant  manufactured  by 
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Figure    5.3       Instrumentation    Plan    of    Test    Structure 


60 


Products  Research  Corporation  and  labelled  1422A1,  the  last 
digit  being  a  1  hour  work  life. 

The  locations  of  the  free  field  pressure  gage  are  shown 
on  the  test  geometry  sketch.  Figure  5.1.   One  is  located  a 
foot  above  the  center  of  the  plate  and  suspended  by  elastic 
shock  cord  and  manila  line.   A  good  view  of  this  arrangement 
is  provided  by  the  photo  in  Figure  5.4.   The  location  of 
this  gage  was  selected  so  that  it  would  sense  the  same 
pressure  wave  which  loads  the  plate.   This  was  considered  an 
important  enough  criteria  to  locate  the  gage  above  the 
plate,  even  though  the  reflected  positive  wave  from  the 
plate  surface  would  reload  the  gage  at  approximately  0.4 
msec  . 

The  second  pressure  gage  should  therefore  be  located 
away  from  the  test  structure,  preferably  out  in  the  free 
field,  on  an  arc  of  the  same  radius  as  the  plate.   Rigging 
limitations  prevented  this  arrangement,  so  it  was  positioned 
at  the  same  radius  as  the  plate  by  clamping  it  to  the  side 
of  the  test  panel.   A  mounting  block,  three  inches  thick  was 
fixed  to  the  panel  so  the  gage  would  be  somewhat  removed 
from  the  boundary  effects  that  the  incident  pressure  wave 
would  create  at  the  fixed  edges  of  the  panel . 

The  pressure  gages  used  were  one  quarter  inch  tourmaline 
crystal  gages,  manufactured  for  the  Navy  at  the  Naval 
Surface  Weapons  Center.   They  are  rated  for  a  reliable 
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Figure  5.4   Test  Rig  Fully  Instrumented 

sensing  of  pressures  up  to  10,000  psi,  although  greater 
pressures  have  been  accurately  measured  with  this  type  of 
gage.   The  Naval  publication  listing  the  specifications  for 
the  tourmaline  gage  [Ref.  15]  indicates  that  the  quarter 
inch  diameter  gage  is  an  adequate  size  gage  for  the  test. 
As  presented  on  page  16  of  the  report,  the  method  of 
determination  of  proper  gage  diameter  is: 

9,  for  explosive  charge  =  .15  msec 
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t„,  transit  time  of  shock  wave  across  gage  =  .25  inch/ 

60,000  in/sec 
=  4.167  xlO"*^  sec 
1.5tp,,  to  account  for  oil  filled  boot  =  6.25  x  10"^  sec 
9/1. 5t^  =  24 

and  from  the  chart  on  page  17, 
R   =  the  pressure  response  ratio 

=  the  apparent  pressure/actual  peak  pressure 
=  .9793 

This  is  well  within  the  acceptable  correction  limit  of  15%. 
The  pressure  gage  greatly  affects  the  accuracy  of  the  peak 
pressure  measurement  due  to  its  finite  size  and,  therefore, 
orientation  of  the  gage  is  important.   The  publication 
indicates  that  all  measurements  and  calibrations  of  the 
tourmaline  gages  are  based  on  a  sideways  or  horizontally 
positioned  gage.   It  was  not  until  after  the  test  was  it 
realized  that  the  second  gage,  located  on  the  side  of  the 
box  in  a  vertical  orientation,  was  not  in  the  ideal 
posi  t ion  . 

After  the  test  structure  was  positioned  in  the  water, 
the  explosive  charge  was  taken  out  to  it  and  the  detonators 
inserted.   The  charge  was  run  down  to  depth  and  secured  in 
position  by  a  lanyard  in  the  rigging.   This  insured  a  center 
position  over  the  plate.   The  eight  pound  TNT  charge  for 
which  all  preshot  calculations  were  made  was  ordered  by  the 
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WCSF.   Although  it  is  considered  a  standard  stock  item  for 
the  Navy,  it  proved  difficult  to  obtain.   Moreover,  when  it 
was  received,  its  physical  condition  was  judged  to  be  so 
deteriorated  that  it  was  destroyed.   A  counter-proposal  was 
made  to  use  what  was  available:   sixteen  1/2  pound  TNT 
charges,  which  could  be  bundled  together  to  form  one  charge. 
This  configuration  is  shown  in  the  lower  right  hand  side  of 
Figure  5.1.   It  was  fixed  with  an  exploding  bridge  wire 
(EBW)  detonator  inserted  in  each  of  the  four  center  charges. 
The  test  was  conducted  on  February  29th,  1984.   Figures  5.5 
through  5.9  show  the  shot  sequence  and  the  two  post  shot 
views  of  the  structure. 

A  listing  of  the  instrumentation  and  the  equipment  used, 
the  recording  set-up  and  a  wiring  schematic  may  be  found  in 
Appendix  E.   All  of  the  equipment  listed  is  the  property  of 
the  West  Coast  Shock  Facility  and  is  maintained  there. 
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Figure  5.5   Lowering  the  Rig  into  the  Bay 
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Figure  5.6   Test  Rig  in  Position 
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Figure  5.7   Water  Plume  from  Shot 
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Figure  5.8   Post  Shot  View  of  Test  Rig 
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Figure  5.9   Post  Shot  Damage 


69 


VI .   RESULTS  AND  INTERPRETATION  OF  DATA 

This  chapter  will  present  the  physical  results  of  the 
test,  then  the  data  obtained,  followed  by  a  discussion  on 
the  limitations  of  the  instrumentation,  then  the  lessons 
learned  from  the  test  and  methods  of  improving  future  tests 
of  this  nature. 

Large  deformations  of  the  plate  were  achieved  and  the 
pattern,  as  seen  in  Figure  6.2  is  as  anticipated: 
symmetrical,  with  the  greatest  deformation  in  the  center 
measuring  1.45  inches.   EPSA  predicted  1.24  inches  as  the 
maximum  displacement  at  the  center  node  at  a  very  early 
time,  1.8  msecs.   A  precise  measurement  of  the  final 
deformation  was  performed  after  the  shot.   A  deformation 
measurement  rig  was  designed  by  the  model  maker,  C.  Crow,  at 
NPS  which  would  support  a  dial  indicator  depth  gage  as  it 
was  traversed  across  the  plate.   The  readings  obtained  are 
presented  in  Figure  6.1  on  the  following  page.   These  were 
taken  before  the  plate  was  removed  from  the  backing 
structure  to  preclude  any  relaxing  that  may  occur  when  the 
plate  is  removed  from  its  rigid  support. 

EPSA ' s  results  are  considered  very  good.   It  must  be 
understood  that  EPSA  did  not  take  into  consideration  any  of 
the  possible  secondary  effects  which  may  have  loaded  the 
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POST  SHOT  MEASUREMENT 
IDENTIFICATION  NUMBERS 
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4)  1.122  15)  1.334 

5)  1.186  15)  0.887 
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Figure    6.1      Post    Shot    Deformation    Measurements 
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Figure  6.2   Panel  Deformation 

plate,  most  notably  the  afterflow  loading  generated  by  an 
accretion  of  water  particles  on  the  plate.   Moreover,  some 
unpredicted  plate  responses  also  occurred  which  otherwise 
would  have  allowed  the  plate  to  deform  to  a  greater  extent. 
There  was  an  unexpected  shearing  of  the  stiffeners  from  the 
edges  at  each  end,  as  seen  in  Figure  6.3.   The  other 
response  that  was  not  planned,  but  considered  a  likely 
occurrence  should  the  loading  be  excessive,  was  a  tearing  of 
the  plate  at  the  center  cavity.   A  close-up  photo  of  the 
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Figure  6.3   Stiffener  Sheared  at  End 

tear  is  shown  in  Figure  6.4.   The  only  other  physical  damage 
to  the  test  structure  was  a  parting  of  two  shackles  which 
attached  wire  ropes  to  one  of  the  pneumatic  fenders,  a 
severe  deformation  of  the  brass  terminal  connection  box  on 
the  side  of  the  backing  structure,  and  a  shearing  of  the 
threads  on  the  hardened  plastic  connector  in  the  base.   Some 
of  these  effects  can  be  seen  in  previous  Figure  5.9. 

Another  view  of  the  deformed  plate.  Figure  6.5  shows  the 
deformation  pattern  and  the  stiffeners  with  only  a  slight 
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Figure  6.4   Tear  Along  Edge  of  Center  Cavity 

out-of-plane  twist.   It  is  evident  that  the  tripping  effect 
desired  did  not  take  place.   Also  to  be  noticed  in  the  photo 
are  that  a  few  of  the  strain  gages  became  detached  from  the 
stiffeners  during  the  test.   One  of  them  is  seen  lifted 
above  the  stiffener  surface.   Not  shown  is  the  shearing  of 
the  lead  wires  which  were  laid  on  top  of  the  stiffeners; 
they  parted  at  the  same  instant  each  of  the  stiffener  ends 
sheared  . 
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Figure  6.5   Stiffener  Deformation 

After  the  shot,  an  accurate  EPSA  prediction  of  stresses 
in  the  plate  was  finally  produced.   Figure  6.6  shows  the  Von 
Mises  stress  levels  for  the  quarter  plate  model  at  time  1.2 
msec  after  the  shock  wave  loads  the  plate.   The  maximum 
stress  is  55305  psi,  well  above  the  yield  stress  of  43200 
psi .   It  is  located  at  the  top  of  the  center  cavity, 
suggesting  that  the  plate  should  tear  first  at  this  point  of 
maximum  stress.   The  experimental  results  confirm  this. 

A  shearing  of  the  stiffener  at  the  ends  is  not  indicated 
by  the  stress  plot.   To  understand  why  requires  a  description 
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Figure  6.6   Von  Mises  Stresses  Predicted  by  EPSA 

of  the  way  EPSA  treats  stiffened  elements.   To  model  the 
stiffener  running  through  a  column  of  elements,  EPSA 
determines  the  mass  of  the  1  inch  by  3/16  inch  stiffener  and 
adds  the  corresponding  amount  to  the  shell  elements  in  the 
column.   The  orthotropic  nature  of  a  stiffened  element  is 
handled  in  the  following  manner.   For  each  of  these 
elements,  strain  increments  are  used  to  compute  stress 
resultants  in  the  shell.   These  same  strain  increments  are 
transformed  to  the  centroid  of  the  stiffener.   Using 
const i tui t ive  equations  for  a  beam,  these  strains  are 
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converted  to  stresses  in  the  stiffener.   Finally,  these  stress 
resultants  are  used  to  modify  the  membrane  forces  and  bending 
moment  in  the  plate  element  in  the  direction  which  the 
stiffener  runs.   The  final  result  is  a  single  stress  at  the 
midthickness  of  the  plate.   With  the  current  method  of  EPSA 
output,  stresses  cannot  be  evaluated  at  the  extreme  fibers  of 
a  stiffener,  nor  at  the  fixed  boundary  where  shearing 
occurred.   It  should  be  noted,  however,  that  Weidlinger 
Associates  have  developed  EPSA  II  which  does  output  the 
necessary  stresses  throughout  a  stiffener.   It  also  allows  for 
the  modeling  of  several  different  stiffener  types. 

Success  was  achieved  in  that  data  was  recorded  on  all 
twenty-four  channels.   As  to  the  quality  and  the  worth  of  the 
data,  that  is  still  to  be  determined.   A  description  of  the 
appearance  of  some  typical  recordings  are  contained  herein, 
along  with  reproductions  of  the  traces,  but  there  is 
insufficient  room  in  this  paper  to  include  all  data  histories. 

Data  reduction  at  the  test  site  consisted  of  a  tape 
playback  into  a  Visicorder  unit.   All  of  the  strain  gage 
histories,  pressure  and  accelerometer  histories  were 
displayed  on  visicorder  paper,  along  with  the  calibration 
signal  for  each  of  the  above  and  a  one  kHz  time  signal. 
Appendix  E  describes  the  instrumentation  and  recording 
procedure.   Tape  playback  was  at  7  1/2  ips  and  Visicorder 
speed  was  80  ips.   This  was  for  an  initial  look  at  the  data 
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and  for  later  reference.   Xerox  copies  were  made  of  the 
Visicorder  sheets  because  they  fade  when  exposed  to  light. 
Next,  the  shock  facility  instrumentation  technician  brought 
down  to  NPS  one  of  the  Ampex  1300  tape  recorders  and  the 
magnetic  tapes.   The  tapes  were  played  back  into  an  HP-5451C 
Fourier  Analyzer  and  recorded  onto  a  hard  disk.   This 
allowed  for  a  local  analyzing  of  the  data,  graphing  on  the 
HP  plotter  unit/  and  fitting  of  the  response  curves  with  the 
time  code  signal  and  calibration  levels.   Most  importantly, 
transfer  of  the  data  to  the  Fourier  Analyzer  permitted  a 
selective  processing  in  the  frequency  domain  and  the 
permanent  storage  of  the  records  at  NPS. 

The  recorded  histories,  from  the  two  pressure  gages 
appear  as  though  the  ratings  of  the  instrumentation  may  have 
been  exceeded.   Also,  it  appears  as  if  there  is  a  signal 
driving  the  instrumentation  beyond  its  maximum  setting  of 
10,000  psi  and  causing  its  output  to  maintain  this  peak 
throughout  the  time  of  shock  wave  loading.   Figure  6.7  is  a 
reproduction  of  the  Visicorder  output.   Shown  are  the  traces 
for  pressure  gages  one  and  two,  located  one  foot  above  the 
panel  and  on  the  side  of  the  panel  respectively,  and  for  the 
1  kHz  time  signal.   Calibration  of  the  amplitude  in  the  time 
domain  was  performed  by  a  simple  procedure  of  laying 
measurements  on  a  3  in.  x  5  in.  card  from  the  calibration 
signals  and  the  1  kHz  time  signals,  then  transferring  them 
to  the  pressure  history. 
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Figure    6.7      Visicorder    Output 
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At  the  test  site,  the  recording  charge  amplifiers  were 
set  to  a  maximum  of  10,000  psi--exactly  what  seems  to  be 
recorded  on  both  gages.   As  shown  earlier,  theory  predicted 
an  incident  pressure  from  an  eight  pound  spherical  TNT 
charge  to  be  3950  psi .   One  would  suspect  that  the  actual 
pressure  developed  by  the  use  of  16  one-half  pound  charges 
may  have  been  somewhat  greater.   Post  shot  inspection 
revealed  physical  damage  to  both  pressure  gages.   The  leads 
inside  the  oil-filled  sensing  boot  had  become  detached  on 
one  gage  and  the  four  tourmaline  crystals  in  the  other  gage 
were  delaminated.   The  time  at  which  this  damage  occurred  is 
unknown  and  such  damage  would  have  caused  erratic  readings. 

Other  damages  for  the  excessive  amplitudes  include  the 
possibility  of  a  calibration  error  in  the  instrumentation 
or  ,  as  reported  in  the  Compendum  of  Underwater  Explosions 
Research  [Ref .  15] ,  the  multiple  charges  have  a  much 
enhanced  effect  over  a  normal  ,  homogenous  charge  of  equal 
weight.   The  pressures  from  each  of  the  charges  may  be 
additive  when  they  are  detonated  sympathetically.   There  is 
much  research  left  to  be  done  in  this  area  of  what  is  called 
the  multiple  charge  effect.   For  example,  the  pressure 
generated  by  a  single  one-half  pound  charge  and  sensed  at 
eight  feet  away  is: 

P  =  22505  (0.5-'-/^/8)  ^*  ^^  =  1473  psi 
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When  using  sixteen  of  these  charges  and  they  are 
arranged  in  such  a  manner  that  sympathetic  detonation  is 
induced,  16  x  1473  =  23,569  psi.   This  is  enough  pressure  to 
saturate  amplifiers  and  damage  pressure  gages.   In  view  of 
this,  some  of  the  questions  raised  include  whether  or  not 
some  or  all  the  charges  were  detonated  sympathetically,  or 
if  there  could  have  been  a  delayed  effect  due  to  the 
geometry  of  the  TNT  block.   The  author  fully  believes  that 
the  shape  and  height  of  the  plume  was  very  different  from 
that  of  a  standard  spherical  charge.   There  may  have  been  a 
concentrating  or  jet  effect  of  the  pressure  due  to  the 
stacking  and  the  geometry  of  the  charges.   To  quote  directly 
from  Cole,  page  229,  "...it  can  be  expected  that  charges 
that  do  not  have  spherical  symmetry  will  give  rise  to  a 
shock  wave  which  is  not  symmetrical,  and  differences  in  form 
of  the  wave  at  different  points  around  the  charge  are  in 
fact  observed ." 

This  statement  leads  one  to  question  what  shape  of 
pressure  history  curve  would  be  generated  from  the 
rectangular  block  of  TNT  used  during  the  test.   Fluid  non- 
linearities  also  affect  the  propagation  of  a  shock  wave  and 
are  difficult  to  predict,  model,  or  reproduce.   As  an  added 
note,  the  equivalent  weight  of  TNT  required  to  produce 
10,000  psi  incident  pressure  at  eight  feet  is  65  pounds. 
Questions  such  as  those  raised  above  can  only  be  answered 
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through  further  expar imentat ion  and  proper  pressure  field 
measurements.   Cole  concludes  his  chapter  on  pressure  waves 
by  stating  that,  "while  an  exponential  curve  is  a  simple  and 
convenient  approximation  to  the  form  of  an  underwater  shock 
wave,  it  is  by  no  means  a  perfect  representation,  and  in 
some  circumstances,  is  a  rather  poor  one."   The  question  of 
what  was  the  true  pressure  loading  experienced  by  the  plate 
may  never  be  known.   The  data  records  are  presented  in  the 
report  for  possible  future  insight  and  understanding. 

Without  the  knowledge  of  the  loading  history  of  the 
plate  an  accurate  prediction/comparison  of  the  strains 
induced  in  the  plate  is  impossible.   A  best  fit  correlation 
was  attempted  between  the  recorded  strains  and  those 
predicted  by  EPSA  for  the  two  possible  loading  extremes;  the 
smallest  being  an  eight  pound  charge  and  the  greatest  being 
a  sixty-five  pound  charge.   The  next  pages  show  these 
strain  history  comparisons  for  several  gages  using  an  eight 
pound  charge.   The  details  of  how  these  strain  histories 
were  generated  by  EPSA  are  included  in  Appendix  C. 

As  can  be  seen  from  the  strain  comparison  (Figures  6.8 
through  6.11) ,  nothing  generated  by  the  computer  seems  to 
come  close  to  agreeing  with  the  recorded  data.   All  twenty- 
one  of  the  strain  records  are  difficult  to  interpret.   There 
is  an  excessive  amount  of  high  frequency  noise  imposed  on 
what  may  be  a  valid  strain  signal  within  the  first  few 
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Gage  8 
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Note:   scales  not  same 
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Figure  6.8   Gage  8  -  Recorded  Strain  and  EPSA  Prediction: 


Figure  6.9   Gage  11  -  Recorded  Strain  and  EPSA  Predictions 
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Figure  6.10   Gage  16  -  Recorded  Strain  and  EPSA  Predictions 
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Fiqure  6.11   Gage  19  -  Recorded  Strain  and  EPSA  Predictions 
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milliseconds.   However,  this  noise  could  not  be  filtered  out 
electronically  by  the  HP-5451C  Fourier  Analyzer.   Curve 
fitting  could  be  attempted,  but  to  even  do  this  would  lead 
to  suspect  results.   Erroneous  signals  were  surely  generated 
as  the  lead  wires  to  the  gages  were  pulled  and  twisted. 
They  were  glued  to  the  plate  surface  which  was  undergoing 
large  deformations  in  very  early  times.   Also,  the  lead 
wires  which  ran  across  the  stiffeners  were  parted  when  the 
stiffeners  sheared.   Moreover,  severe  damage  was  done  to  the 
leadwire  connection  box  on  the  side  of  the  backing 
structure;  wires  were  stripped  and  broken  as  the  shock  wave 
opened  the  box  to  the  environment.   Recommendations  to  avoid 
some  of  these  problems  in  future  tests  are  included  later  in 
this  chapter.   Appendix  F  contains  select  strain  histories. 
Presented  in  Table  II  are  the  results  of  the  strain  gage 
analysis,  with  no  conclusions  drawn.   Only  the  amplitudes  of 
the  first  two  peaks  are  tabulated;  the  validity  of  the 
subsequent  recording  is  questionable.   It  is  thought  that 
the  first  sharp  peak  in  the  strain  histories  is  the  firing 
pulse  signal.   This  is  not  shown  in  Figures  6.8  though  6.11, 
but  may  be  seen  in  the  data  located  in  Appendix  F. 

The  one  instrument  record  which  is  readable  and  shows 
good  results  is  the  accelerometer  history.   This  is  shown  in 
Figure  6.12  with  the  integrated  velocity  and  displacement 
records.   The  HP-5451C  Fourier  Analyzer  was  used  to  perform 
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TABLE  II 
Strain  Readings  from  Gages 


Firing 

Rise 

First 

Gage 

Pulse 

Time 

Duration 

Peak 

Number 

(  ustrain) 

(msec) 

msec) 

(  ystrain) 

1 

-2250 

0.09 

0.11 

+  2000 

2 

-1440 

0.07 

0.10 

+  1200 

3 

-1650 

0.10 

0.12 

+  750 

4 

-400 

0.05 

0.10 

+  800 

5 

-1350 

0.07 

0.10 

+  950 

6 

-750 

0.06 

0.  11 

+  1200 

7 

-2800 

0.08 

0.10 

+  1300 

8 

-800 

0.08 

0.  12 

+  1400 

9 

+  500 

0.05 

0.05 

-1000 

10 

+  1200 

0.04 

0.06 

-2100 

11 

+  2200 

0.06 

0.10 

-2500 

12 

+  800 

0.05 

0.07 

-2200 

13 

+  800 

0.03 

0.06 

+  1200 

14 

* 

* 

* 

* 

15 

-1200 

0.07 

0.09 

+  1000 

16 

-12000 

0.04 

0.07 

+  2000 

17 

-1200 

0.03 

0.04 

-1000 

18 

-4750 

0.07 

0.10 

-1000 

19 

-3000 

0.04 

0.08 

+  700 

20 

* 

* 

* 

-1000 

21 

* 

* 

* 

* 

Amplitude  too  small  to  read. 


No  attempt  was  made  to  preserve  the  sign  (+  or 
during  transcipt ion  . 
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Figure  6.12   Accelerometer  Histories 
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the  integrations.  The  three  records  appear  as  they  should: 
the  integration  of  an  impulse  is  a  step  and  the  integral  of 
a  step  is  a  ramp  function. 

Each  of  the  instrument  records  in  Appendix  F  clearly 
shows  the  firing  pulse  impressed  on  the  electronic  signal. 
This  indicates  that  the  wiring  configuration  was  such  that 
it  was  particularly  susceptible  to  noise  pickup  and 
electrical  interference.   The  first  singular  voltage  spike 
is  identified  as  the  firing  signal  because  it  appears  in 
each  of  the  records  in  similar  form.   It  is  followed  by  a 
"dead-time"  on  each  record  where  no  other  signals  occur  for 
at  least  1.8  msec.   This  time  corresponds  to  that  necessary 
for  the  shock  wave  to  travel  nine  feet  from  the  charge  to 
the  plate.   The  shape  of  the  firing  pulses,  and  of  many  of 
the  other  sharp  peaks  on  the  records,  are  too  similar  and 
cause  one  to  question  why.   An  explanation  is  found  in  the 
Shock  and  Vibration  Handbook  by  Harris  and  Crede  [Ref .  16] . 
It  presents  a  thorough  discussion  on  how  improper  or 
inadequate  sensing,  recording,  and  playback  equipment  will 
limit  the  frequency  of  signals  which  can  be  accurately 
measured.   For  example,  listed  on  page  17-15  of  the 
reference,  the  upper  limit  of  FM  magnetic  tapes  is  5000  Hz, 
a  light  beam  galvanometer  is  limited  to  50,000  Hz,  an  AC 
carrier  amp  with  demodulator  is  5000  Hz,  and  the  same  for  an 
AC  power  supply.   These  numbers  are  representative  only,  and 
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the  equipments  are  not  the  same  as  those  used  during  the 
test.   However,  the  characteristic  rise  time  for  an  impulse 
shock  wave  is  five  or  six  milliseconds  (as  limited  by  the 
tourmaline  1/4  inch  gage) .   This  means  that  the  electronic 
equipment  must  be  sensitive  to  200  kHz.   The  tape  recorders 
used  were  Ampex  1300 's  which  are  rated  at  20  kHz  as  a 
response  time.   This  is  precisely  the  recorded  rise  times 
for  the  firing  pulses  and  peak  signals  on  the  data  records. 
Table  II  also  lists  the  rise  times  for  these  signals.   It  is 
evident  that  the  limitations  of  the  equipment  were  exceeded 
during  this  shot.   What  should  be  used  for  tests  of  this 
nature  are  FM  tape  recorders  with  a  much  better  response 
time,  such  as  an  Ampex  2200  or  a  Honeywell  101.   A  note  of 
caution,  however,  is  that  the  higher  the  frequency  of  the 
recording,  the  poorer  the  signal  to  noise  ratio. 
Considerations  must  be  given  to  a  proper  arrangement  of 
filters  and  playback  speeds  to  obtain  the  best  records. 
Oscilloscope  photography  has  provided  good  response 
histories  in  past  experiments.   This  method  should  always  be 
considered  a  possibility. 

To  properly  set  up  the  electronics  for  an  underwater 
explosion  test  is  an  art  in  itself.   Many  lessons  were 
learned  and  some  of  them  will  be  briefly  mentioned  as  items 
that  should  be  considered  in  any  future  tests: 

-There  must  be  a  zero  time  signal  on  one  of  the  data 
channels.   This  is  necessary  during  the  analysis  phase  so 
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the  time  code,  the  gage  histories,  and  the  zero  time  can  be 
lined  up  to  measure  small  differences  in  arrival  times. 

-If  clipping  of  the  signals  is  present,  ensure  that  it 
is  not  caused  by  the  playback  electronics  being 
oversensitive . 

-For  the  tourmaline  pressure  gage  calibration,  the 
Q-step  method  described  by  Cole  (pages  182-3)  is  highly 
recommended . 

-There  are  corrective  techniques  to  back  out  actual 
peaks  in  signal  records  should  they  be  rounded  due  to 
equipment  limitations.   One  of  these  methods  is  discussed  by 
Ronald  Tussing  in  his  text  on  page  33  [Ref.  14]. 

-The  firing  circuit  should  be  coupled  to  a  transformer 
so  as  not  to  be  grounded  with  the  other  gages.   In  this 
test,  magnet ic/ induct ive  pickup  through  the  cables  is 
thought  to  have  caused  the  firing  signal  to  be  impressed  on 
the  other  lines  and  perhaps  creating  a  signal  of  such 
magnitude  as  to  overdrive  the  amplifiers.   In  future  tests 
it  is  essential  that  the  cable  used  to  transmit  the  firing 
voltage  be  positioned  as  far  away  from  the  recording  lines 
as  possible . 

A  chapter  in  Harris  and  Crede's  text  is  titled, 
"Interference  and  Noise  in  Transmission  Cables".   It  brings 
out  many  important  considerations  and  the  means  to  correct 
them.   For  example,  on  page  12-21,  it  defines  electrical 
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interference  or  pickup  as,  "Noise  components  superimposed  on 
the  desired  signal  due  to  the  proximity  of  the  connecting 
cable  to  the  electromagnetic  field  of  an  electrical 
disturbance"...  "and  ' static- type '  disturbances  from 
switching  transients  are  the  worst  offenders."   A  solution 
offered  is  the  use  of  properly  shielded  cables  and  the 
correct  grounding  of  the  circuit  in  relation  to  the 
electronic  equipment. 

Another  source  of  signal  interference  mentioned  by  the 
text  comes  from  the  movement  of  the  signal  carrying-cables. 
On  page  12-26  it  cautions  that,  "Noise  is  generated  when  the 
cable  is  suddenly  squeezed,  bent,  struck,  or  mechanically 
distorted.   Peak  noise  voltages  from  this  source  were 
frequently  as  large  as  the  actual  acceleration  signals  being 
recorded."   A  shock  test  is  a  dynamic  phenomena  and  lead 
wires  to  instrumentation  will  be  twisted  and  pulled  as  the 
shock  wave  impinges  on  them  and  as  the  plate  deforms.   There 
seems  to  be  almost  no  way  to  eliminate  this  problem,  however 
there  may  be  methods  to  reduce  the  effect  this  has  on  the 
signal.   Rather  than  fixing  the  strain  gage  lead  wires 
directly  to  the  plate  with  the  PRC  adhesi ve/ sealant ,  it  is 
recommended  that  a  sheath  or  tube  arrangement  be  devised  to 
permit  the  lead  wires  to  move  as  the  plate  deforms.   This 
should  be  done  at  least  in  the  immediate  area  of  the  test 
panel  where  large  displacements  will  occur. 
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Additional  considerations  must  be  given  to  the  use  of 
the  correct  type  of  strain  gage.   There  are  special  high- 
elongation  gages  which  can  measure  up  to  10%  strain  with 
accuracy.   Also  pointed  out  in  Reference  IG  is  that  the 
dynamic  type  gage  will  have  the  greater  gage  factor  and 
provide  maximum  possible  electrical  signal  for  a  given 
strain.   The  use  of  an  iso-elastic  foil  will  have  a  greater 
gage  factor  than  constantan  and  therefore  a  greater 
resistive  charge  for  a  given  strain. 

Another  important  aspect  of  the  use  of  strain  gages  is 
that  they  are  bonded  properly  to  the  test  specimen.   As 
pointed  out  by  Harris  and  Crede  on  page  17-5,  the  proper 
functioning  of  a  strain  gage  is  completely  dependent  on  the 
bond.   One  must  use  the  manufacturer's  recommended  cement 
and  follow  the  directions  closely.   Incorrect  strains  will 
result  if  the  bond  is  not  completely  over  the  entire  area  of 
the  gage,  or  if  it  becomes  partly  or  fully  detached  during 
the  test  . 

The  recording  of  transient  strains  under  dynamic  loads 
is  a  mostly  difficult  task,  as  brought  out  by  the  work  at 
M.I.T.  by  Dr.  Witmer  [Ref.  17].   In  a  series  of  shock  tests 
performed  at  the  Aeroelastics  and  Structures  Research  Lab, 
only  about  0.2  msec  of  data  was  the  best  that  could  be 
obtained.   Although  the  severity  of  the  impulse  for  the 
experiments  was  much  greater  than  that  for  the  underwater 
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shock  test,  the  environment  was  much  more  controlled.   The 
problems  most  often  encountered  were  gage  detachment  and 
lead  wire  rupture.   A  very  helpful  procedure  for  the  bonding 
of  strain  gages  is  presented  by  Dr.  Witmer  in  his  report, 
particularly  with  regard  to  the  use  of  backing  material  to 
cover  the  gage.   His  procedure  is  included  in  Appendix  D  for 
future  reference. 

As  a  final  caution  in  this  discussion  on  the  instrumen- 
tation of  a  shock  test,  it  would  be  wise  to  thoroughly  check 
out  the  electronic  recording  set-up  just  prior  to  the  test 
so  that  it  is  well  known  what  will  be  the  charactec  of  the 
signals,  as  distinct  from  interference  noise.   Suggested  is 
a  simple  procedure  of  turning  off  all  power  supplies  and 
striking  the  plate  to  see  the  strain  gage  response  within 
the  elastic  regime.   It  is  also  most  important  for  the 
individual  conducting  the  test  to  have  a  complete 
understanding  of  both  the  nature  of  the  test  (including 
underwater  shock  phenomena,  material  response  of  the  test 
structure,  the  desired  results,  etc.)   and  the  electronic 
measuring  and  recording  of  the  data.   Both  of  these  go  hand- 
in-hand  and  must  be  tailored  to  provide  meaningful  data. 

One  final  recommendation  for  future  work  in  this  area  is 
a  redesign  of  the  test  panel  to  allow  for  better  tripping  of 
the  stiffeners.   A  single  stiffener  in  the  longer  direction 
of  the  test  section  will  provide  a  greater  area  for  a 
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stiffener  to  suffer  compression.   Also,  investigation  into 
the  effects  of  detaching  the  ends  of  the  stiffener  so  it  is 
free  to  fall  flat  would  be  of  interest.   In  any  case,   it's 
hoped  that  the  work  to  date  has  paved  the  way  for  further 
research  in  this  field  of  study. 
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APPENDIX  A 
UNIAXIAL  TENSION  TEST  DATA  FOR  6061-T6  ALUMINUM 

Uniaxial  tension  tests  were  performed  to  characterize 

the  material  properties  of  the  6061-T6  Aluminum.   A 

rectangular  section  of  the  material  was  cut  from  the  end  of 

the  blank  and  six  round  test  specimens  were  machined  from 

it.   Figure  A.l  shows  the  dimensions  of  a  specimen.   The 

specification  used  was  ASTM  E8-69  for  small-sized  round 

tension  test  specimens  proportioned  to  a  standard  2.0  inch 

gage  length.   The  charts  from  the  six  tests  are  reproduced 

as  Figures  A. 2  through  A. 7.   Of  the  six,  only  four  tests  are 

considered  valid.   On  test  #1,  the  chart  speed  was  too  slow; 

it  was  increased  from  2.0  inches  per  minute  to  20.0  while 

the  drive  speed  was  kept  at  0.2  inches  per  minute.   For  test 

#2,  the  original  diameter  was  not  measured  accurately  and 

the  results  are  only  approximate.   Tests  #3  and  #6  are  good 

and  the  results  are  listed  below. 

Outside 
Test      Diameter 
Number    ( inches) 

1 

2  est.  0.250 

3  0.251 

4  0.255 

5  0.251 

6  0.248 

Average 

of  last  four         43200  44900  11.50 
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Yield 

Ultimate 

Elongation 

Stress 

Yield  Strength 

from  Chart 

(psi) 

(psi) 

(%) 

44600 

46300 

12.3 

43250 

44866 

10.68 

42980 

44742 

12.04 

43148 

44866 

12.04 

43267 

45026 

11.36 

O 

In 


O 

t-' 

<r 

\xl    Ul 

CC   UJ 

<^  X 

o 

o  3 

q: 

Lu     O 

S 

^ 

H-    - 

^    1 

o: 

v/l      J> 

^ 

o^      _- 

o 

a: 

\r, 

vj 

O 

O 

a: 

4J 

UO 

Z 

O 

U- 

o 
o 

,-*! 

"o 

o 
m 


o 


CO 

c 

0 

•H 

cn 

c 

(U 

£ 

■H 

Q 

C 

^ 

0) 

\o 

e 

o 

■H 

o_ 

0 

4.   t 

Si 

___    ^   1 

a 

'o 

en 

yr\ 

r< 

■p 

X 

' 

en 

."^N 

\ 

(D 

^          ^ 

. 

Eh 

\.__/, 

1 

y 

0 

- 

•H 

M 

en 

o 

c 

o 

(U 

Eh 

—  +  1 

v=. 

<o 

.H 

m 

oo    ^ 

* 

< 

^--^    — ^ 

lr^ 

0 

u 

0 
IT" 
•H 

fo 

98 


-p 
en 

0) 


+J 

U 
n3 

u 

c^ 
c 

•H 

o 
u 
d) 


<: 


99 


tim-^ 


M 
(1) 

Eh 


u 
(d 

u 
c 

•H 
73 
U 
O 
U 
(U 


ro 


•H 


100 


1 

T-- 

-  - 

-1 

,p- 

T--- 

: 

, 

1 

..:i: 

i 
-1  -  - 

!■■■ 

;i   .. 

*'> 

11  3 

! 

.■  1 

:t- 

• !  ■ 

.  .1.." 

:; 

\ 

! 

1 

-L-- 

-.  1 , 

°  ■     '■s 

•i_.  .  -  — 

— 

-  — 

::  ..I-, 

:5]^ 

V  - 

1 

:i' 

'i 

i  '  f 

, 

V- 

-, — 

--—+ 

-— 

■y- 

: ' : .  i ; :  ■ ' 

:l 

I-.-: 

'  ■ 

■  ~" 

~"^~' 

-  — 

i 

— 

-  - — 

-~-r 

— 

f 

- — 



— 

— 

r— 

-— - 

'■•-■ 

!':• 

-— 

mI 

'■:■ 

— : 

::  1 

""/ 

i 
[ 

— 

L_ 
i 

:':. 

■' 

-— 

1 

— 

j  ■ 

;' 

— 

iM 

i  '■  ■  ■ 

.,!-. 

— 

-• 

— 

irrr: 

\      : 

— 

'  1  ^ 

-- 

-- 

'"'. 

— 

- 

— 

.... 

i 

;' 

f  ■■■ 

'['I 

!— 

U- 

111. 

j 

■  ■  r ; 

' ' ; . 

.'■■ 

'■A 

■..''. 

:iLhL":| 

j! 

—ff 

1 

i'i' 

j : ' : 

!;' 

v3 

: '.  1 

:^'i 

•  i  i  ■ 

1; 

1 : 

:~ll 

" 

"]~\ 

;:'' 

!  ■  r 

'  ■ ! ! 

*'* ', 

■  .■  1 

!'  i  ■ 

Cb 

■* 

\'"' 

.  :  ,  j 

i     j 

■  1 

' '  1 

iil^ 

'  .  .Ul 

ill' 

1  ;  ' 

~" 

'— ■ 

1  . 

— - 

— ;• 

!  1 ' 

-^ 

,  ^ 

f 

^— 

l'; 

— 

--— 

■ii' 

'  !  1  ' 

M    i 
1 

\\U 

' '  ■ . 

i':: 

I ' '  ■ 

:  r!  i 

'  ■  I 

.: 

'.;i 

";;: 

1  - 1 ,. 

■V 

,1" 

t:.- 

:'- 

-y\. 

1 

"  i 

::.: 

:-" 

'^\ 

■  i ' ' 

1 !  r' 

;    :! 

,.'! 

i : 

:  ■-   i  ■   : 

'  -  '  i 

:   1 "  ■ :  ■ 

1  ._-, 
1 

j 

1 

■ii 

.    .    L  L 

■   ! 

X 

-  — 

— 

. 

1 ,  -.  I 

■    i    ■    ' 

1 

I  ■  • : 

::.. 

— - 

■-- 

.-- T 

F': 

i  j 

■ 

— 

:F' 

:  ■■:[ 
■Ii! 

:|: 

:■■: 

::., 

1 

--'7 

i  ■  ■ '  ' ; . ; 

1 

:  '  ■ 

'  i  ;i 

;     ■  1 

!  1  !  i 

: 

1  '  ' 

•    1  - 

f  '  ■'  ( 

o 

, 

' 

'"h": 

I  ; 

■  ■  \ 

''; 

';■' 

— 

1     T 



— 

— - 

— 

. — 

":"r~ 

■11 

J  !   ,   ' 

: 

^'1 

:; 

:!' 

1          ,   ; 

r4^ 

' ;  i  ; 

:  '  1  1 



!:_.! 

_._. 

■  J  ' 

.:- 

r^i 

--- 

._ 

.._ 

— 

.... 

■I;l;:': 

.. 

■          '■'■■'\ 

\  \ 

.  ' 

'  J, 

H 

—  1   -  - . 

--3— 

i 

:  ;  i  : 

■:  t 

b 

■- 

--- 

-- 

u:. 

':  ■ 

'' 

:   :■   ■::  | 

•  1 
'-■■4 

■  .  !  ■ 

■'i! 

::,;!.  1_ 

:.:  ' 

!|- 

n 

■■  1 

' 

i: 

.::|,: 

■1 

^v: 

... 

■  i  i 

■-':' 

■  1   M 

, 

. 

, 

, 

, 

\| 

":',l 

•-:■ 

-:- 

■'/ 

— ,  :  -,- 

-:: 

"i    \ 

'-::- 

S 

-  - 

--- 

— 

.... 

'ii: 

{ 

■  '  : 

"kit 

1 

■  ( ■ 

1 

^ 

^ 

i        :i:| 

ll'.A 

!-- 

__  . 

— - 

-— 

1 

-_ 

^ 

^' 

__:. 

.  I  (  : 

;  , 

, 

;  I  i 

u-L' 

::[: 

i 

"i 

" '  ■ 

.' 

i 

nX 

v^ 

1 

;: 

-■ 

*n 

- 

; 

' 

'^ 

^N 

i 

— - 

■a 

■ ' ; 

;:.- 

1 

-^- 

-;--- 

-  - 

.... 

~ 

-- 

- 

-- 

V 

: "  i ' 

■"     1 

</) 

!   '\ 

1 

■   1 

■  ■?= 

:       1 

CJ 

1 

1 

1 

1 

1 

.■  1 

-P 


-P 

u 

u 

en 

c 

•H 
T! 
U 
0 
O 
(U 


< 
U 

cn 

■H 


101 


(U 

Eh 


u 

u 

c^ 
c 

0 

o 


in 

< 

en 


102 


in 
w 

Eh 


td 
u 

D^ 
C 
•H 

-r) 

0 
o 

0) 


IsO 


u 

o 

•H 

Cm 


103 


-p 
en 
CD 


-P 

S-i 

to 

x: 
u 

c 

•H 

13 

O 

u 

CD 


QJ 
U 

a 

•H 


104 


APPENDIX  B 
DESIGN  DRAWINGS  OF  THE  BACKING  STRUCTURE  AND  TEST  PANEL 

The  drawings  used  for  the  manufacture  and  assembly  of 
the  test  structure  at  the  Naval  Postgraduate  School  are 
found  on  the  following  pages.   Also  shown  are  the  sketches 
used  to  order  the  plate  cut  to  size. 
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Figure  B.l   Test  Panel  and  Backing  Structure 
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APPENDIX  C 
EPSA  ON  THE  VAX 

The  EPSA  User's  Manual  titled  "Underwater  Shock  Response 
of  Submerged  Shells"  [Ref.  9]  presents  the  cotnmands 
necessary  to  execute  an  EPSA  job  on  a  VAX  computer. 
Francois  Daube's  thesis  [Ref.  10]  also  shows  how  to  run  EPSA 
at  NPS  and  how  to  generate  and  view  the  PATRAN  files.   The 
intent  of  this  appendix  is  to  cover  some  of  this  material 
again  and,  hopefully,  bring  out  some  helpful  information 
that  was  previously  omitted. 

The  input  deck  used  for  the  pre-shot  analysis  is  shown 
again  on  the  following  page  as  Figure  C.l,  the  name  of  this 
input  data  file  is  CPLT.DAT.   Most  of  the  input  is 
adequately  described  in  Reference  9  ;  those  modifications  to 
EPSA  which  deal  will  the  PATRAN  interface  are  not  found  in 
the  reference,  but  are  described  in  Figure  C.2. 

To  access  PATRAN,  two  numbers  at  the  end  of  line  2  must 
be  added  to  the  input  deck.   This  specifies  the  number  of 
PATRAN  files  to  be  created.   The  only  limitations  on  the 
numbers  are  that  they  be  positive  integers  and  do  not  exceed 
the  number  of  time  steps  requested.   However,  it  is  advised 
to  keep  them  small  due  to  the  additional  memory  allocation 
these  files  require. 
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NPS  PLATE-  2  REC  STIF   STRN  TOP 
4000   1000   100   0044 
.0000005     1   4000   100 
1   1   17   1   16 

16  10   4   1   1   0   27   20   26   14   1   1   0 

0  0 

10000000.       .3       .000255      40000.       .18750      0,0      0.00001 
'stif    10      1 
-111 
187      0 

1  0      0      0.0      0.0 

17  1      0      9.0      0.0 
-17    17      10      0.0      0.6 
0/ 

160      1      1      16      10 

11      17      01010 

2  1      11      10010 

3  1      17      11110 

4  1      11      11110 

1         1      1      1      0.0      0.0      0.0 
0.0      0.0      0.0      0.0      0.0      0.0      0.0      0.0/ 
111112113      17    11      17    12 
17    13671672673      17    11    1 
17    11    2      17    11    3      12    7       1      12    7      2      12    7      3 

I  11  11  11  21  11  3611612 
613  12  11  12  12  12  13  17  51 
17    5      2      17    5      3/ 

111112114115  16  11 
16  12  16  14  16  15671672 
6  7  4  6  7  5  16  10  1  16  10  2  16  10  4 
16    10    5      11    10    1      11    10    2      11    10    4      11    10    5/ 

II  5   10  10  1   9   11  1   10  12  1   9   12  10  10 
16  5   9   16  5   10  16  1   9   16  1   10  11  1   10 

10  1  10  12  8  9  11  8  10  10  8  9  16  1  7 
16  1   8   11  1   11  11  2   11  11  3   11  11  4   11 

11  5   11  11  6   11  11  7   11  11  8   11  11  9   11 
11  10  11/ 

113      17    13      17    11    31       11    3573 

5  11  3   12  2   3   12  7   3   12  11  3   17  7   3 
11  1   3   11  3   3   11  5   3   11  7   3/ 
0.5625    0.1875   1.000 

11       1/ 

0/ 

9.0    0.0000945    60000.0 

1    108.0    0.0     9.0 

8.0    22505.0     1.18    0.053    -0.135   0.0008/ 

0    0    0    0    0 

Figure    C.l      EPSA    Input    Deck    for    Pre-Shot    Analysis 
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liad  of  line  2:             A  4 

number  of  PaCran  number  of  Pacran 

displacement  plots  stress/strain  plots 

End  of  line  3:            1  AOOO             100 

timestep  number  time  step  number    number  of 

to  scare  time  to  stop  plotting    increments 
history  plot 

Entire  line  4:           1  1           17      1             16 

number  of  nodes  number  of      node  identifer     element 

to  be  plotted  elements      first  row, number  17   numtjer 


Figure  C.2   Modifications  to  the  Input  Deck 

The  numbers  added  to  the  end  of  line  three  will  cause 
EPSA  to  create  a  file  of  data  points  during  a  run.   The 
subroutines  PLOTEN  and  BIGPLOT  were  added  to  EPSA  to  do 
this.   These  data  points  may  be  plotted  on  the  Tektronics 
4013  terminal  by  using  a  separate  Fortran  program   to 
produce  either  a  displacement  or  strain  history.   If  this  is 
requested  by  adding  the  numbers  to  line  three,  then  an 
entire  line  four  must  also  be  input.   The  data  files  created 
are:   FOR022.DAT  for  displacements,  FOR021.DAT  for  Y 
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strains,  FOR024.DAT  for  X  strains,  and  FOR023.DAT  for 
velocities.   Although  four  files  are  created,  a  program  to 
access  the  velocity  history  at  a  node  (the  FOR023.DAT  file) 
has  not  been  written.   To  access  the  displacement  and  strain 
history  files,  the  following  Fortran  programs  were  written: 
PL0TD2,  PLOTSX,  and  PLOTSY.   Listings  of  these  files  are  at 
the  end  of  this  appendix.   To  plot  one  of  these  histories 
from  a  Tektronics  terminal,  simply  type  the  command,  "RUN 
PL0TD2"  (or  PLOTSX  or  PLOTSY)....  and  return. 

An  example  of  a  strain  history  is  presented  in  Figure 
C.3,  and  the  data  file  which  created  it  is  found  in  Figure 
C.4.   However,  it  is  the  executable  version  of  these  files 
which  is  run.   These  are  created  by  compiling  the  Fortran 
file,  then  linking  to  PLOT  10  by  the  following  commands: 
$  FOR  PL0TD2.F0R 
$  LINK  PL0TD2.0BJ  'PLOTIO' 

Any  changes  made  to  these  routines  must  be  supported  by 
the  PLOTIO  software  package--see  a  PLOTIO  User's  Manual  for 
details . 


115 


t-\ 
o 

u 

•H 

c 
o 

+J 

M 
Q) 

Eh 

c 

O 

0) 

+J 

o 

a, 

>i 

o 
■p 

CO 

•H 

a: 


(T3 


m 
U 

S-i 


116 


.s  Q  \ 


p£&< 


■,-i^ 


r^'i" 


-  ^ 

.05JQ3J5C.-,  , 

~  -4 

.  3  »  0  5  u  S  3  E  -  -  5 

I 

.1^  =  ?!:J1E-0  4 

1 

.  0  4  T  u  1  =  ■,  E  -  1  4 

1 

.4aso^7ij.,ia 

-•5 

.o?77,?-4E-0u 

-1 

.  35770? '^E-n-^ 

-1 

.577Q033S-03 

.3 

.5?0?S21E-03 

-2 

.4817341E-03 

.  3 

.2336Pbr£-05 

-P 

.  o  P  ?  2  7  ?  7  E  -  n  5 

-5 

.=■61  375T£-.'!5 

-5 

.  :  1 ?02a4£.05 

•  -  5 

.3?'4  05'7E-"3 

-5 

.0?4a, 04E-OJ 

-4 

.3^61 3o5E-03 

-4 

.4*3  =  ^41  IE- 'I? 

-'4 

.2?«5444E-0 J 

-  5 

.  :?77i4iE-D3 

.  J 

.J01P350E-'i5 

-1 

. qqajosir -0  5 

-1 

.■^55JS°2E-0  5 

-  J 

.44) ojaaf-A^ 

-? 

.  4o56o4''E-03 

-4 

.  4  0  0  ■»  2  "  3  c  -  0  5 

-5 

.3342b03E-'i3 

-3 

.52733°lE-03 

.  ;a 

.2304?'oE-0  3 

-1 

,  lua; 335F-03 

-1 

.3e567^Sc-0? 

.  j 

S1-j05!4C-np 

"^ 

^  4  7  ;  i  4  7  ;  .  ,;  3  > 

-1 

00-i?[i '.'-'? 

-  1 

5a?'^i3  =  E-0? 

-1 

ooi 0 1 1 nt- :? 

-1 

l„aap3=c-02 

.a  _ 

55t3oq,£-03 

-=(  . 

07335=; ''E-n 

-5  . 

550<373aE-.1j 

-a. 

31 4674o£-05 

-a. 

535=730£-03 

-5  . 

300a54?E-03 

.  3  _ 

°e?b5oi£-0J 

-1  . 

sSfca3SuE-0? 

-3  . 

3Ami,0  3E-OP 

.J  . 

7  3  0  5  4  1  :  E  -  ■!  P 

-  I 

;C7a:;,  ac.,lP 

'^-  >. 

SM~7^---.i  P^ 

-  5  . 

3j^ai  ji.;.i)p 

-  1  . 

-  1  3  "  ?  1  1  £  -  ■"!  ? 

S  . 

-'?75  3  7lc-A-i; 

1  . 

a^37;7PF--P 

■-2j 

4  0  u  5  M  T '  - 1  r^ 

?. 

J  a  p  q  3  0  0  £  -  0  P 

1 . 

£a9067P£-OP 

^. 

4  0  0a0  30E-'"i3 

I . 

l0  =  ''446E-i'>? 

1  . 

l^S  =  •)4,•l£-.1P 

5^- 


•t-)-'  - 


I  .35 
I  .22 
1.11 
1  .IS 
=  .23 
=  .53 
3.50 
1.01 
1  .2'J 

1  .54 


1.33 
1  .35 

?.:i 

?.  10 


324-. 
4023 

4  73  I 

5=12 
71  45 
1  2'^^ 
552- 
3°t  1 

PU45 

1  7aj 
52aS 
4873 

43o0 
aoQ  1 

41  5o 
S:ia2 

5773 
923" 

d7TP 

5  =  35 
1  020 

4  =  35 

7045 

57sa 
0051 
aa5'' 
a  a  a  a 
5«20 
al62 
5S03 

00,4 

05  7  1 
0725 
3147 
533- 


2E--I3 
3E-'i5 
3E-03 
=  £-■13 
2  £  -  •"'  4 
7£-'lJ 
7E-03 
2  £  -  0  3 

1  £  -  0  5 
JE-03 

2  £  -  0  4 

1  £  -  ">  4 

0  r  .  '1  J 

-!  E  -  0  4 

1  £  -  0  u 
1  E  -  0  3 
2E-''>4 
3£-0a 

3  £  -  0  3 
3E-0  3 
4E-0  5 
3E-'-'3 
OE-OP 
7E-02 
3  F  -  0  ; 
5E-02 
4E-02 
9  E  -  0  3 
?E  • 
7£- 


3E  - 

3£  - 
3£  - 
3C  . 


0  2 

OP 

0  2 
sE-  12 

J  F  -  0  p 


Figure    C.4       Strain    History    Data    File,    FOR024.DAT 


117 


After  the  creation  of  the  input  deck,  the  comnnands 
required  to  execute  the  EPSA  job  are  shown  in  the  figures 
below.   Figure  C.5  shows  the  queries  and  responses  necessary 
for  an  interactive  EPSA  run.   Note  that  the  request  for  a 
Virtual  Mass  Array  filename  is  returned  anpty.   This  is  done 
only  for  the  case  of  a  flat  plate  model  where  only  the  Plane 
Wave  Approximation  is  used.   If  a  cylinder  were  being 
modeled,  the  proper  response  would  have  been  EPM.VMA. 
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Another  comment  about  the  interactive  runs  is  to  be  aware 
that  interactive  runs  should  be  made  only  if  the  user 
intends  on  executing  only  a  small  number  of  time  steps  and 
wants  to  remain  at  the  terminal  during  the  run.   Otherwise, 
a  batch  job  is  more  appropriate.   This  also  has  the 
advantage  of  allowing  the  user  to  access  and  work  with  his 
other  files  during  the  run;  the  terminal  screen  is  not  busy 
wiwth  the  EPSA  run.   The  commands  for  the  batch  execution 
are  shown  in  Figure  C.6. 
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In  Figure  C.7,  note  that  the  time  history  and  the  core 
tape  file  names  are  requested.   The  file  names  input  are 
consistent  with  those  names  used  by  the  authors  of  EPSA. 
The  time  history  file  would  be  used  for  the  plotting  of 
output  data.   The  specific  plotting  routine  to  do  this  is 
not  available  at  NPS.   An  example  of  its  use  may  be  found  in 
Reference  11.   The  core  tape  file  is  used  for  a  restart  of 
EPSA--if  it  was  desired  to  start  the  computation  process  at 
the  same  timestep  which  a  job  terminated.   An  example  of 
this  procedure  is  found  in  the  User's  Manual.   Another 
feature  of  the  batch  process  is  illustrated  in  Figure  C.7. 
A  file  name  BEPSAI.LOG  is  generated  during  execution,  but  it 
cannot  be  accessed  until  the  job  is  terminated.   Should  a 
request  to  view  (type)  this  file  be  made  before  job 
termination,  the  DCL  symbol  $  is  returned,  meaning  "not 
ready  yet".   This  is  a  positive  signal  that  the  job 
submission  has  gone  through  and  computation  is  in  progress. 
The  BEPSAI.LOG  file  contains  only  system  control  messages 
from  the  VMS  level ,  such  as  the  date  and  time  completed  and 
CPU  time  required  for  execution. 

An  example  of  other  DCL  commands  not  covered  adequately 
in  the  VAX/VMS  literature  is  shown  in  Figure  C.8  on  the 
following  page.   A  request  is  made  to  display  all  the  files 
in  the  user's  directory  of  name  FOR018.DAT  with  the  version 
unspecified  (;*) .   A  second  example  follows  the  first.   Note 
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that  the  EPSA  run  had  created  the  previously  unmentioned 
files  FOR015.DAT  and  FOR017.DAT.   These  two  data  files  are 
not  used,  they  are  in  machine  language,  and  their  creation 
was  meant  to  be  eliminated  by  Francois  Daube . 
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To  access  PATRAN,  follow  the  procedure  in  Figure  C.9. 
The  work  TEK  is  input  when  using  a  Tektronics  graphics 
terminal  and  it  must  be  in  upper  case  letters.   The  screen 
will  clear  and  follow  with  what  is  shown  in  Figure  CIO. 
Answer  with  GO.   The  other  commands  and  options  necessary  to 
utilize  the PATRAN  graphics  are  adequately  covered  in 
Francois  Daube's  thesis.   Refer  to  reference  15.   However, 
one  feature  of  the  PATRAN  package  not  discussed  is  illus- 
trated in  Figure  C.ll.   To  stop  the  plotting  of  a  view  which 
is  not  desired  or  incorrect,  simply  hit  the  Control  and  C 
keys  simultaneously.   This  will  abort  the  function  without 
losing  the  PATRAN  session  files.   A  control  Y  will  cause 
an  exit  from  PATRAN  altogether  and  a  loss  of  all  session 
files  . 

Should  it  be  necessary  to  make  changes  in  EPSA,  or  any 
other  executable  file  on  the  VAX,  the  procedure  is  shown  on 
page  29  of  the  User's  Guide  [Ref.  16].   An  example  is  now 
given.   The  required  changes  would  be  made  in  the  edit  mode 
to  the  Fortran  version  of  the  file  or  subroutine.   An 
example  would  be  to  change  the  line  CX  =  -CX  to  read  CX  =  CX 
in  the  subroutine  CPTSTRN.   Now  the  corrected  version  of 
EPSA5804.FOR  must  be  compiled.   To  do  this  without  a  request 
for  debugging  type 

$FORTRAN  EPSA5804.FOR 
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An  object  file  is  created  named  EPSA5804 . OB J .   For  ease 
in  making  future  changes  to  EPSA,  a  library  file  named 
EPSA5804.OLB  should  be  created.   This  is  done  by  the  command 

$LIBRARY/CREATE  EPSA5804.0LB  EPSA5804.OBJ 

Next  the  library  file,  or  the  object  file  if  none  was 
created,  must  be  linked  to  the  main  program  to  produce  an 
executable  file.   Type  in 

$LINK  EXE=EPSA5804I.EXE  EPSAI , EPSA5804/LIBR 

if  linking  to  the  library,  or 

$LINK  EXE=EPSA5804I .EXE  EPSAI , EPSA5804 . OB J 

if  none  was  created.   The  purpose  and  value  of  a  library 
file  is  that  if  one  were  created  and  another  change  be 
necessary  to  subroutine  CPTSTRN,  the  entire  EPSA  code  would 
not  have  to  be  recompiled  and  linked.   Simply  "cut  and 
paste"  the  subroutine  from  EPSA5804.FOR  so  that  it  is  a 
separate  Fortran  file,  CPTSTRN.FOR.   Make  the  necessary 
corrections,  then  compile  this  smaller  file  to  produce 
CPTSTRN.OBJ.   Now  link  this  object  file  to  the  EPSA  library 
and  it  will  replace  the  older  version  of  the  routine.   A 
form  of  the  command  which  achieves  this  is 

$LINK  EXE=EPSA5804I .EXE  EPSAI ,CPTSTRN , EPSA5804/LIBR 
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Of  course  there  will  not  be  a  correct  Fortran  file  of 
the  current  executable  version  (EPSA580_I . FOR)  within  the 
directory.   Either  a  printout  must  be  kept  updated  through 
manual  changes,  or  the  Fortran  file  on  the  computer  should 
be  changed  but  not  compiled. 

If  a  change  were  made  to  one  of  the  Fortran  files  which 
plot  the  output  histories  on  the  Tektronics  terminals,  the 
compiling  and  linking  is  achieved  by  typing 

$FOR  PL0TD2 

$LINK  PL0TD2, ' PLOTIO • 

PlotlO  is  the  plotting  package  that  contains  commands 
used  by  these  routines.   As  mentioned  earlier  in  this 
chapter,  to  plot  one  of  these  histories  on  the  Tektronics, 
type  in  RUN  PL0TD2. 


131 


APPENDIX  D 
STRAIN  GAGE  INSTALLATION  TECHNIQUE 

The  procedure  used  for  this  test  is  presented  here  as  a 

reference  for  future  work  and  possible  improvement.   The 

procedure  proved  successful  as  all  gages  indicated  some 

strain.   It  is  based  upon  the  BLH-SR-4  application  kit  (part 

103158)  manufactured  by  BLH  Electronics,  42  Fourth  Avenue, 

Waltham,  MA   02254.   The  figures  used  are  from  the  BLH 

instruction  sheet.  The  materials  required  are 

-Permabond  910  catalyst  P/N  215203  in  jar  with  brush 
applicator  (red) 

-Permabond  910  P/N  223012  in  squeeze  bottle  (clear) 

-Sheet  of  teflon  film 

-Tweezers,  razor  blades,  silicon  carbide  400  grit 
sandpaper 

-Scotch  Magic  Brand  transparent  tape 

-Tissues  or  Kimwipes 

-Rubbing  alcohol  and  acetone  in  plastic  squirt  bottles 

A.   SURFACE  PREPARATION 

1)  Clean  general  gage  area  to  be  bonded  so  it's  free  of 
dirt  and  grease . 

2)  Scribe  lightly  crosshairs  over  the  center  position, 
using  care  not  to  score  severely  or  it  will  create  a  stress 
concentrator.   Use  a  lead  pencil  if  you  don't  want  to  use  a 
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metal  scribe,  but  remember  that  it  will  be  erased  during  the 
cleaning  process  and  must  be  redone  before  the  laying  of  the 
gages . 

3)  Use  silicon  carbide  paper  to  smooth  out  the  surface. 
Wipe  clean  with  Kimwipes. 

4)  Saturate  a  tissue  with  acetone  and  rub  well.   Now, 
use  alcohol  to  remove  the  acetone. 

5)  Repeat  the  above  steps  until  you  are  certain  that 
the  surface  is  absolutely  clean  and  free  of  dust  and 
oxidation.   The  wiping  tissue  should  not  be  discolored. 

B.   BONDING  INSTRUCTIONS 

1)  Using  clean  tweezers,  remove  a  strain  gage  from  its 
package  and  lay  it  down  in  the  exact  position  where  it  is  to 
be  bonded.   Do  not  touch  the  back  of  the  gage  with  your 
fingers  or  allow  it  to  be  contaminated  in  any  way.   Most 
strain  gages  have  crosshairs  embedded  in  the  backing 
material  which  are  specifically  made  for  this  alignment 
process.   For  future  use,  cut  a  section  of  teflon  film  about 
2  inches  by  3  inches  and  have  it  available. 

2)  Remove  a  length  of  scotch  tape  about  3  inches  and 
carefully  fix  one  end  of  it  down  on  the  surface  and 
straighten  it  out  so  that  the  tape  will  extend  beyond  the 
gage  about  one  inch.   Now  lay  the  tape  down  evenly 
proceeding  from  the  fixed  end  until  the  gage  is  now  picked 
up  by  and  securely  stuck  to  the  tape.   Double  the  tape  back 
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Figure     2 


Figurt     3 


Figure     4 


Figure     5 


Figure     6 


Figure     7 


Figure  D.l   Strain  Gage  Bonding  Technique 
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over  itself  to  form  a  loop  and  secure  the  other  end  of  the 
loop  to  the  plate.   This  tape  will  be  used  to  replace  the 
gage  in  the  exact  position  after  the  adhesive  is  applied. 

3)  Apply  a  catalyst  to  the  back  of  the  gage  sparingly, 
yet  ensure  the  entire  gage  is  covered.   It  should  dry 
quickly . 

4)  Now  lift  up  the  working  end  of  the  tape  and  re-lay 
it  so  the  gage  is  over  its  final  position,  but  not  quite 
touching  the  surface.   While  holding  the  tape  with  one  hand, 
position  the  precut  piece  of  teflon  film  so  that  it  is  now 
between  the  back  of  the  tape  and  the  fingers  you  will  use  to 
rub  the  tape  after  the  adhesive  is  applied  (this  is  the  same 
hand  as  the  one  still  holding  the  tape  above  the  surface) . 
This  can  be  tricky. 

5)  Apply  the  catalyst  by  placing  only  a  few  drops  onto 
the  surface  close  to  where  the  tape  is  now  attached  so  that 
it  will  spread  through  the  area  where  the  gage  will  lay. 

Be  careful  not  to  touch  the  tip  of  the  applicator  to 
anything  or  it  will  contaminate  the  entire  contents. 

6)  Now  continue  laying  down  the  tape  so  that  the  gage 
contacts  the  surface,  first  only  one  end  of  it,  then 
proceeding  to  the  other  end.   Permabond  will  ooze  out  from 
under  the  surface.   This  is  the  reason  for  the  Teflon  film-- 
to  prevent  your  fingers  from  becoming  glued  together.   Use 
the  film  to  rub  the  backing  of  the  tape  to  make  a  good  bond 
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between  the  gage  and  the  surface,  and  press  down  on  the  gage 
for  45  seconds. 

7)  Allow  the  perinabond  to  set  about  5  minutes,  then 
carefully  peel  up  the  working  end  of  the  Scotch  tape.   Pull 
backwards  and  to  the  side  at  a  ninety  degree  angle  to  the 
line  of  the  tape.   Work  the  tape  back  until  it  is  entirely 
removed.   Be  particularly  careful  that  you  do  not  exert 
enough  force  on  the  edge  of  the  gage  to  cause  it  to  detach 
in  the  process . 

8)  Now  the  gage  is  bonded  securely  in  place.   Keep  in 
mind  that  during  this  procedure  the  accuracy  of  the  strain 
gage  is  entirely  dependent  on  a  complete  and  thorough 
bonding  to  the  surface  which  is  to  be  measured. 

9)  After  attachment  of  the  lead  wires  to  the  gage 
terminals,  an  application  of  "Gage  Cote"  to  the  surface  of 
the  gage  will  provide  additional  protection.   This  is 
advisable  in  a  workshop  environment  where  the  test  plate  is 
exposed  to  rough   treatment. 

C.   IMPROVED  METHODS 

In  the  study  of  impulsively  loaded  structures  performed 
by  E.A.  Witmer  at  MIT,  an  improved  method  of  strain  gage 
application  was  devised  and  reported  as  Appendix  A  of 
reference  2.   Although  the  test  arrangement  was  quite 
different  from  that  studied  in  this  paper,  the  gage 
installation  technique  proved  successful  for  the  measurement 
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of  transient  strains.   Dr.  Witmer  presents  an  "enhanced 
survival"  method  of  protecting  the  gages  under  severe 
impulse  loading  conditions.   To  prevent  the  gages  from 
becoming  detached  from  the  surface,  a  covering  patch  of 
polyimide  sheet  was  cemented  over  the  entire  gage  and  a 
portion  of  the  lead  wires.   Figure  D.2  is  a  reproduction  of 
his  illustration,  followed  by  his  description  as  found  on 
page  101  of  reference  2. 

"From  a  0.001  inch  thick  sheet  of  polyimide,  cut  a 
'cover  patch'  of  the  size  indicated  in  Figure  A-7  (Figure 
D-2  here)  .   Roughen  and  clean  this  patch  as  in  Step  5  (using 
a  400  grit  silicon-carbide  paper  and  clean  with  N5) .   Apply 
a  thin  coat  of  AE-15  cement  to  (a)  the  gage  area,  (b)  the 
polyimide  insulator  patch,  and  (c)  both  sides  of  the 
polyimide  cover  patch.   Position  as  shown  in  Figure  A-7  and 
cover  the  entire  region  with  a  strip  of  0.003  x  1  x  1  inch 
of  teflon  film.   Secure  the  teflon  in  place  with  cellophane 
tape.   Repeat  Step  8  to  complete  the  cementing/curing 
process  (this  step  describes  the  curing  of  the  bonded  gage 
with  proper  pressure  and  heat) .   The  teflon  cover  does  not 
adhere  and  is  hence  removed.   The  result  is  a  multilayer 
cover  of  cement/polyimide/cement ,  with  an  exterior  layer  of 
AE-15  cement." 

Dr.  Witmer  used  special,  high  elongation  strain  gages 
manufactured  by  Micro-Measurement  Co.,  Romulus,  Michigan. 
His  gage  data  follows: 
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GAGE  EP-08-125AD-120 
(.40"L  X  .22"W) 


POLYIMIDE    COVER   PATCH 
(7/8"   X   5/8"    X    .001") 


POLYIMIDE    INSULATOR   PATCH 
(3/8-   X   3/8"    X    .001") 


GAGE   EP-08-031DE-120 
(.12-    X    .27") 

POLYIMIDE    COVER  PATCH 
(5/8"    X   1/2"    X    .001") 

POLYIMIDE    INSULATOR   PATCH 
(3/8"   X   3/8"    X    .001") 

FIG.    A. 7       POLYIMIDE    COVER   PATCH    DIMENSIONS    AND    LOCATION 


Figure  D.2   Improved  Strain  Gage  Mounting  [Ref.  2] 


Gage 

length  (in.) 
without  tabs 

length  (in.) 
with  tabs 

width  (  in  . ) 

Polyimide  backing 
length  ( in . ) 
width  ( in . ) 
thickness  (in.) 


EP-08-031DE-120 
0.031 

0.140 

0.032 

0.27 
0.12 
0.001 


EP-08-125AD-120 
0.125 

0.250 

0.125 

0.40 
0.22 
0.001 
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special  consideration  to  lead  wire  arrangement  is  also 
given  in  reference  2.   In  an  effort  to  make  the  lead  wires 
as  nearly  massless  as  possible,  very  fine  (0.003  inch 
diameter)  varnish  covered  copper  lead  wires  (AWG  40)  were 
used.   They  were  suspended  above  the  structure  in  a  "soft 
coil"  arrangement  to  reduce  the  forces  on  the  wires  and  the 
connection . 

Although  a  configuration  such  as  this  appears  necessary 
to  properly  measure  transient  strains,  this  particular 
method  is  not  feasible  for  external  gages  in  an  underwater 
shock  test.   Attention  must  be  given  to  this  matter  in 
future  tests.   As  recommended  earlier,  a  tube  or  sheath 
arrangement  to  permit  the  lead  wires  to  extend  seems 
plausible  . 
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APPENDIX  E 
INSTRUMENTATION  AND  EQUIPMENT  AT  WCSF 

The  following  list  identifies  the  electronic  equipment 
used  in  the  instrumentation.   The  components  listed  are 
located  at  the  West  Coast  Shock  Facility,  San  Francisco,  CA. 

Unit  Quantity 

Ampex  FR  1300  tape  recorder  2 

Ampex  magnetic  oxide  coated  tape  2 

General  Radio  Co.  1952  Universal  Filters 
(high  pass/low  pass) 

Bell  and  Howell  DC  amplifer  CEC  1-168  16 

Endevco  signal  conditioner  model  4470  15 

Current  Regulated  Bridge  Conditioners 

Systron  Donner  Time  Code  Reader/Search  Unit  8134        1 

Honeywell  Universal  Bridge  Balance  Unit  2 

Model  #B2-6,  six  channels  each  unit 

Honeywell  1858  Visicorder  1 

Firing  Circuit,  800  volts  DC  1 

Relay  operated,  capacitive  discharge 

Dupont  EBW  Detonators  X-175  4 

Tourmaline  pressure  gages,  1/4  inch,  four  pile  2 

with  a  50  foot  cable  attached 

Micro  Measurements  120  OHM  strain  gages  21 

EA-06-250BG-120 

Endevco  piezores i stive  accelerometer  1 

model  2262-2000,  range  +2000  g 
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Seacon  Connector,  24  pin  XSM-CCP  with 
50  foot  cable,  Boston  insulated  wire 
and  cable  co .  MWF-19  NEOP 

Seacon  Penetrator ,  XSM-BCR 


Bellfoil  Shielded  Cables, 


Belden  8424, 
Belden  9776, 


1/4  inch 
1/2  inch 
16  gage 


1-200  feet 
2-200  feet 


Standard  U.S.  Navy  Brass  Connection  Boxes 
one  on  side  of  backing  structure 
one  floating  50  feet  from  the  test  rig  and 
200  feet  from  pier 


141 


LftKiD'-iNT   C/^B'^'i 


"voicr 


FiRiT    liC-JaMMST-S 


5  R 1  Otf.  r 


AmPLiF' t~RS 


J_ 


■DC 


^Tir^E"  <^^0£ 


r 


I 


Ar^PViF'iTRS 


VliiCoRbt^ 


Figure  E,l   Flow  Chart  for  Electronic  Setup 


142 


APPENDIX  F 
SELECTED  DATA  RECORDS 
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